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ABSTRACT 

We have studied multi-band surface brightness profiles of a representative sample of 34 nearby dwarf 
irregular galaxies. Our data include Galaxy Evolution Explorer (GALEX) FUV/NUV, UBV, Ha 
and Spitzer 3.6 ^m images. These galaxies constitute the majority of the LITTLE THINGS survejl^. 
By modeling the azimuthal averages of the spectral energy distributions with a complete library of star 
formation histories, we derived the stellar mass surface density distributions and the star formation 
rate averaged over three different timescales: the recent 0.1 Gyr, 1 Gyr and a Hubble time. 

We find that, for 80% (27 galaxies) of our sample galaxies, radial profiles (at least in the outer 
part) at shorter wavelengths, corresponding to younger stellar populations, have shorter disk scale 
lengths than those at longer wavelengths, corresponding to older stellar populations. This indicates 
that the star-forming disk has been shrinking. In addition, the radial distributions of the stellar mass 
surface density are well described as piece-wise exponential profiles, and ~ 80% of the galaxies have 
steeper mass profiles in the outer disk than in the inner region. The steep radial decline of star 
formation rate in the outer parts compared to that in the inner disks gives a natural explanation for 
the down-bending stellar mass surface density profiles. Within the inner disks, our sample galaxies on 
average have constant ratios of recent star formation rate to stellar mass with radius. Nevertheless, 
~ 35% (12 galaxies, among which 7 have baryonic mass < 10® M©) of the sample exhibit negative 
slopes across the observed disk, which is in contrast with the so-called "inside-out" disk growth scenario 
suggested for luminous spiral galaxies. The tendency of star formation to become concentrated toward 
the inner disks in low mass dwarf irregular galaxies is interpreted as a result of their susceptibility to 
environmental effects and regulation through stellar feedback. 

Subject headings: galaxies: dwarf - galaxies: irregular ~ galaxies: evolution - galaxies: stellar content 



1. INTRODUCTION 

^ ' Dwarf irregular (dirr) galaxies numerically dominate 
" actively star-forming galaxies in the local universe (van 
\0 : den Bergh 1966; Gallagher & Hunter 1984; Dale et al. 
■ 2009). Compared to local spiral galaxies, dIrr galax- 
] ies are generally smaller, bluer, less luminous, and gas- 
, rich. dIrr galaxies tend to have low metallicity and a flat 
^-H ' metallicity gradient within measurement uncertainties in 
7— I I the interstellar medium (ISM; e.g. Kobulnicky & Skill- 
T— I , man 1996; Hunter & Hoffmann 1999; van Zee & Haynes 
L; 2006; Croxall et al. 2009). 
. !^ The star formation (SF) process in seemingly sim- 
: pie dIrr galaxies is unclear. The classical Toomre grav- 
5_] ' itational instability model (Toomre 1964; Quirk 1972), 
d which seems to explain the SF threshold successfully 
in spiral galaxies (Kennicutt 1989; Martin & Kennicutt 
2001) and elliptical galaxies (e.g. Vader & Vigroux 1991), 
predicts that the gas surface density of dIrr galaxies, both 
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at intermediate radii and outer radii, is usually well below 
the threshold of large-scale SF for a single-fluid gas disk 
(Hunter, Elmegreen & Baker 1998, hereafter H98; Leroy 
et al. 2008). Even though recent studies (e.g. Yang, Gru- 
endl & Chu 2007) suggest that incorporation of the grav- 
itational contribution of the stars may make the gas disk 
marginally stable, the uncertainties (e.g. the gas/stellar 
velocity dispersion, the disk scale height) included in the 
stability analysis are substantial. Nevertheless, widely 
distributed SF, although apparently inefhcient, has been 
observed in dIrr galaxies, which implies that some other 
internal or local process must play an important role (e.g. 
Elmegreen & Hunter 2006). 

H98 found that, for most dIrr galaxies, the current 
SF rate (SFR, traced by Ha) correlates with the optical 
(e.g. V) starlight profiles better than with the usually 
mass-dominant atomic gas surface density. The correla- 
tion with older stars indicates that SF may be a self- 
regulating process and/or the SF may be affected by 
the existing stars to a certain degree. By analyzing opti- 
cal and near-IR (NIR) broadband data. Bell & de Jong 
(2000) and Bell et al. (2000) found that the SF histo- 
ries (SFHs, described as luminosity-weighted ages there) 
primarily correlate with the if-band surface brightness 
in both spiral galaxies and low surface brightness (LSB) 
galaxies. Recently, Leroy et al. (2008) found that the 
SF efficiency (SEE ^ SFR/Matomic gas) shows a strong 
relationship with the stellar mass surface density, rather 
than with the gas, in atomic gas-dominated regions, such 
as dIrr galaxies and the outer disks of spiral galaxies. Shi 
et al. (2011) proposed an "extended Schmidt law" with 
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explicit dependence of the SFE on the steUar mass sur- 
face density. Therefore, an important question to be ad- 
dressed is how do the SF, SFH, and steUar mass surface 
density relate to each other? 

Generally speaking, the SFR, SFH, and stellar mass 
density could be derived directly from analysis of re- 
solved stellar color-magnitude diagrams (CMD). How- 
ever, besides the problem of stellar crowding in in- 
tense star-forming regions, currently full-scale — in terms 
of both the observational depth and spatial coverage — 
CMD analysis of dirr galaxies is still rare (e.g. Cole et al. 
2007; Hidalgo et al. 2011). Obtaining deep CMDs which 
reach the old main seqiience (MS) turnoffs is crucial to 
reconstructing the whole SFH of a galaxy. As pointed out 
by Gallart, Zoccali, & Aparicio (2005), beyond ~2 Gyr 
ago, shallow (e.g. completeness limit M/ < 2) CMDs at 
best can reliably constrain the integrated SF. Therefore, 
analyzing spatially integrated multi-wavelength broad- 
band data is still an efficient and helpful way to constrain 
the SFHs for a large sample of dIrr galaxies. 

Specifically, ultraviolet (UV) emission is mainly sen- 
sitive to the recent (< 100 Myr, e.g. Kennicutt 1998) 
SF; the optical broadlDand emission has long been shown 
to be able to constrain the SF during the past ~ Gyr 
(e.g. Gallagher, Hunter, & Tutukov 1984); the NIR bands 
have the least sensitivity to recent SF and dust extinc- 
tion compared to shorter wavelengths, and thus give the 
tightest constraint on the stellar mass (e.g. Elmegreen 
& Elmegreen 1984; Kendall et al. 2008). Even more, the 
combination of optical and NIR data is able to largely 
break the notorious age-metallicity degeneracy (Worthey 
1994; Anders & Fritze-v. Alvensleben 2003; Lee et al. 
2007). CardicI et al. (2003) found that the inclusion 
of an NIR band improves the constraints on the stel- 
lar populations by ~ 30 times over using optical colors 
alone. Bell & de Jong (2000) found that there are signifi- 
cant variations in stellar mass-to-light ratios even in NIR 
bands, which is consistent with recent modeling results 
(Dutton 2009). Above all, by fitting the whole spectral 
energy distribution (SED) from the UV to the NIR, the 
stellar mass distribution and recent (< Gyr) SFHs can 
be reasonably constrained. 

With the advent of GALEX (Martin et al. 2005) and 
Spitzer (Werner et al. 2004), acquiring multi-band data 
from the far-UV (FUV) to the NIR for large samples 
of dIrr galaxies has become possible. The Spitzer In- 
frared Array Camera (IRAC, Fazio et al. 2004) 3.6 /im 
passband is primarily stellar emission from the mass- 
dominant older populations (e.g. Helou et al. 2004), even 
though there could be a significant contribution from hot 
dust emission in intense star-forming regions (e.g. Zhang, 
Gao, & Kong 2010). Furthermore, in low-metallicity dIrr 
galaxies, the contribution of polycyclic aromatic hydro- 
carbon (PAH) dust emission to the 3.6 /im passband is 
substantially weaker compared to that in high-metallicity 
spiral galaxies (Engelbracht et al. 2008), which suggests 
that the IRAC 3.6 /im band is a reliable tracer of under- 
lying older stellar populations in dIrr galaxies. 

In an attempt to probe the relationship between the 
stellar mass surface density, SFH, and recent SF, this 
work derives stellar mass surface density profiles and 
rough SFH variations as a function of radius by syn- 
thesizing multi-wavelength data from the GALEX FUV 
to Spitzer 3.6 /xm of a representative sample of 34 dIrr 



galaxies. 

This paper is organized as follows: In § 2, we describe 
the sample and introduce the data used in this work. 
The multi-band surface brightness profiles are presented 
in § 3. Then the stellar population synthesis method used 
in this work follows in § 4. In § 5, we present the results 
of the population modeling, including the stellar mass 
surface density profiles and the SFH variations with ra- 
dius. The discussion and the summary are provided in 
§ 6 and § 7, respectively. We also test the reliability of 
our SED modeling in Appendix A. 

2. DATA AND PHOTOMETRY 
2.1. Data 

The 34 galaxies included in this work (listed in Ta- 
ble 1) are drawn from the LITTLE THINGS sample, 
which is a representative sample of 41 nearby dIrr galax- 
ies spanning a large range in galactic parameters: in- 
tegrated luminosity {Mb of —9 to —19), surface SFR 
density (0 - 1.3 M0yr~^kpc^^), (atomic) gas richness 
(0.02 - 5 Mq/Lb), and central surface brightness /io^ 
(18.5 - 25.5 mag arcsec"^). The sample galaxies were se- 
lected to be relatively isolated. LITTLE THINGS is a 
large VLA project that was granted nearly 376 hours of 
VLA time in the B, C, and D array configurations to 
obtain deep H i-line emission maps of dIrr galaxies with 
high angular and velocity resolutions. The subsample of 
34 galaxies studied here was mostly chosen because these 
galaxies have the full complement of data available, in- 
cluding GALEX FUV/NUV. UBV, Ha and Spitzer 3.6 
fim. The Two Micron All Sky Survey (2MASS) JHKs 
images are too shallow to be used in this work. However, 
UGC 8508, which cannot be observed by GALEX due 
to a nearby bright star, and DDO 165 which does not 
have FUV observations, are also included in the current 
sample. Another three galaxies without an Ha detection 
(LGS 3, M81dwA, DDO 210) were also included. 

The acquisition and reduction of the UBV and Ha nar- 
rowband images of our sample were described by Hunter 
& Elmegreen (2004, 2006; hereafter HE04, HE06). The 
GALEX UV data for 21 galaxies included in the current 
sample were presented by Hunter, Elmegreen & Ludka 
(2010). 

There have been significant improvements in the flat 
flelds and absolute calibration of the GALEX imaging 
data from the GR2/3 pipeline to the GR4/5 pipeline. 
The work of Hunter et al. (2010) used the GR2/3 
pipeline data for some galaxies. Therefore, we obtained 
the GR4/5 imaging data of all of our galaxies from the 
GALEX archive. Most of the Spitzer 3.6 /im images were 
obtained and reduced by the Local Volume Legacy survey 
(LVL, Dale et al. 2009) and the Spitzer Infrared Nearby 
Galaxy Survey (SINGS, Kennicutt et al. 2003). Besides 
the standard post-pipeline processing, additional pro- 
cessing, including distortion corrections, rotation of the 
individual frames, bias structure and bias drift correc- 
tions, has been done by the teams of these two surveys. 

2.2. Surface Photometry 

We follow the procedures as laid out in HE04 for the 
surface photometry. First, the background-subtracted 
images were geometrically transformed and aligned with 
the F-band image. Background-and-foreground sources 
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not belonging to the galaxy were masked. Then the sur- 
face photometry was carried out with the Image Re- 
duction and Analysis Facility (IRAPH!]) task ELLIPSE, 
adopting the same geometrical parameters as derived 
by IIE04. Briefly, the center of the galaxy, position an- 
gle, and the ellipticity were determined from a contour 
in the outer half of the V-band image that was block- 
averaged by factors of a few to increase the signal-to- 
noise. The center was fixed as the geometrical center of 
this isophote, and the major axis is the longest bisector 
that passes through the center that, as much as possible, 
symmetrically divides the galaxy. The surface photom- 
etry of the UBV and Ha images has been described in 
HE04 and HE06. 

To perform photometry on the GALEX images, as de- 
scribed in Hunter et al. (2010), we masked foreground 
stars and background galaxies identified mainly in the 
higher resolution optical band images. For the removal 
of sky background, because the background was quite 
flat, we measured the sky level in tens of square regions 
(10 pixel by 10 pixel) sampled along an ellipse around the 
galaxy, but far enough from the galactic center to avoid 
emission from the target galaxy. The final sky level was 
determined as the average of the sky values in all sampled 
subregions. As to the uncertainties of the background re- 
moval, the standard deviation of the mean sky values 
(Jl^sky) among the surrounding sky regions is usually an 
order of magnitude smaller than the mean of the stan- 
dard deviation in individual regions {a sky)- Therefore, 
only the {(Jsky) were considered. 

For the Spitzer 3.6 fiia images, we first applied the 
same photometry mask as was used for the optical band 
images. Background galaxies can significantly affect the 
photometry on 3.6 /im images of faint dwarf galax- 
ies. Whenever available, high resolution archival Hubble 
Space Telescope (HST) imaging data were inspected to 
further remove the background galaxies based on mor- 
phology. 

Finally, by dividing the luminosity differences between 
adjacent ellipses by the area of the annulus, we obtained 
azimuthally-averaged surface brightness profiles. A value 
of 0.15 mag was taken as the calibration uncertainty for 
the GALEX UV data, and 0.1 mag for Spitzer 3.6 /im. 
The calibration uncertainties of the UBV data for each 
galaxy were given in HE06. We adopt a 20% calibration 
uncertainty for all of the Ha photometry. 

Routinely one uses the average Milky Way (MW) ex- 
tinction curve {Rv=Av / Eb-v—'^A, Cardelli, Clayton & 
Mathis 1989) to correct the photometry for the Galac- 
tic extinction. However, there is considerable scatter of 
the extinction curve shape from sightline to sightline 
through the Galaxy (Fitzpatrick 1999), which could sig- 
nificantly infiuence de-reddening of the short-wavelength 
data (e.g. UV, U) when only an average extinction 
curve was adopted. Fitzpatrick (1999) gave the typical 
uncertainties {crE(\-v)/E{B-v)j Table 1 in Fitzpatrick 
(1999)) when correcting the data with the average ex- 
tinction curve. In this work we included the uncertainty 
{E{B - V)xaE(\-v)/E{B-v)) of adopting the average 

tt IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 



MW extinction curve. The final uncertainties assigned 
to each annulus are a quadratic sum of four contri- 
butions: {(Tsky), photometric calibration uncertainties, 
Poisson noise and E{B — V)^^ e(\-v) / e(b-V) - The flux 
calibration uncertainties usually dominate all other con- 
tributions except in the outer disks. 

3. MULTI-BAND SURFACE BRIGHTNESS PROFILES 

The surface brightness profiles of Ha, FUV, B, and 
3.6 /J,m for all galaxies are presented in the first column 
of Figure [T] The multi-band surface brightness profiles 
broadly follow each other. Nevertheless, for about half 
of our galaxies, the surface brightness profiles of Ha and 
FUV peak in the circumnuclear regions or even in the in- 
termediate radii, rather than the center. The recent SF 
indicator Ha (and also FUV) usually follows the shorter 
wavelength stellar emission better than the longer wave- 
length emission, as is quantified below (Table 2; Figure 
[5]). Therefore, the correlation between Ha and the optical 
bands seen by others was partly caused by a dominant 
contribution of recent SF to the optical passbands. The 
three blue compact dwarf (BCD) galaxies (Haro 29, VT 
IZw 403, Mrk 178) included in our sample show obvious 
central enhancements even at 3.6 /im, which suggests 
that the central regions have sustained the central star- 
burst that characterizes BCDs for an extended period of 
time (e.g. Zhao, Gu & Gao 2011). 

To quantify the radial fall-off of the multi-wavelength 
emission, we obtained scale lengths of the exponential 
disk, Rd, of the FUV, NUV, U, B, V and 3.6 /im surface 
brightness profiles by fitting a single or multi-piece ex- 
ponential to the profiles with the least squares method. 
The results are given in Table 2. Negative values of the 
scale lengths in Table 2 indicate that the relevant sur- 
face brightness increases toward larger radii. To visual- 
ize the difference of the scale lengths at different wave- 
lengths, is plotted against for 
the galaxies (Figure [5]). For the double exponential pro- 
files, we only plot the disk scale length measured in the 
outer part. As is shown, the disk scale lengths measured 
at shorter wavelengths are smaller than those measured 
at longer wavelengths for ^ 80% (27) of our sample galax- 
ies, at least in the outer disks. From Table 2. ~ 65% (22) 
of the galaxies have shorter disk scale lengths measured 
at shorter wavelengths in the inner disk, and 6 galax- 
ies with broken surface brightness profiles have shorter 
scale lengths measured at shorter wavelengths in both 
the inner and outer disk. 

FUV traces SF over the past ~ 100 Myr, i3-band emis- 
sion is dominated by stellar populations younger than a 
few Gyr, and 3.6 /im is a very good proxy for SF over 
the whole lifetime of the galaxy (see references in the 
Introduction). The low metallicity of dirr galaxies, com- 
bined with the lack of radial metallicity gradients, in- 
dicates that the multi-band surface brightness profiles 
should be predominantly determined by the distribution 
of the different stellar populations, with negligible effects 
due to changing internal extinction. Therefore, the fact 
that shorter wavelengths have shorter scale lengths most 
probably implies that the star-forming disks are shrink- 
ing towards the inner regions for most of the dIrr galax- 
ies in our sample. Below, we will come back to this point 
with the help of multi-band SED modeling. The trends 
we see here in dIrr galaxies are in striking contrast with 
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the findings in the luminous spiral galaxies (e.g. Ryder 
& Dopita 1994), for which the inside-out disk growth 
scenario has been suggested and the disk scale length 
decreases with age. 

We note that for more than 50% of the sample galaxies 
the multi-band surface brightness profiles show obvious 
breaks, especially at shorter wavelengths. R^l^t) I ^Diiri) 

is plotted against J^j^^^^t)! ^^D(in) ^^"^ galaxies exhibit- 
ing obvious surface brightness breaks in Figure [3] As 
can be seen, most of the galaxies with broken profiles 
are characterized as down-bending profiles. Furthermore, 
for the down-bending profiles, there is a trend that the 
breaks in FUV profiles are stronger (smaller ratio of 
RD(out)l RD{in) for a down-beuding profile) than that in 
3.6 /im, indicating that the recent SFHs may play an im- 
portant role in shaping the broken profiles. We note that 
the smaller ratio of Rr){mt)/^r>(in) '^^ primarily due to 
the much steeper radial decline of the FUV in the outer 
disk for most galaxies. A down-bending profile has been 
found in the majority of spiral (e.g. Pohlen et al. 2002; 
Pohlen & Trujillo 2006) and dirr (HE06) galaxies. 

4. SED MODELING TECHNIQUE 

4.1. Creating a Library of Star Formation Histories 

Detailed studies of high-quality stellar CMDs for very 
nearby dIrr galaxies (e.g. Dolphin 2000; Skillman et al. 
2003; Cole et al. 2007; Tolstoy, HiU & Tosi 2009) sug- 
gest that low-mass dIrr galaxies exhibit a wide variety of 
SFHs over their lifetime. Therefore, rather than adopt- 
ing the commonly used two-component form of SFHs 
for larger galaxies (an exponentially declining compo- 
nent superimposed by random bursts, e.g. Kauffmann 
et al. 2003; da Cunha, Chariot & Elbaz 2008), we use 
multi-component population models to create our SFH 
library. Specifically, we divided logarithmically the life- 
time from the present to 13.70 Gyr ago into six inde- 
pendent age bins (0 - 0.03 Gyr, 0.03 - 0.10 Gyr, 0.10 - 
0.35 Gyr, 0.35 - 1.18 Gyr, 1.18 - 4.02 Gyr, 4.02 - 13.70 
Gyr), where each age bin has a constant SFR. The way 
we choose these age bins reflects the fact that the sepa- 
ration between isochrones of different ages strongly de- 
creases with age. Our choice of the first age bin (0 - 0.03 
Gyr) is justified by the nearly constant UV colors in the 
first ^ 0.03 Gyr of evolution (e.g. Leitherer et al. 1999). 
Similarly, the last age bin (from ~ 4 Gyr to 13.7 Gyr) 
is justified by the nearly un-evolving shape of the opti- 
cal to NIR spectrum at ages older than ~ 4 Gyr (e.g. 
Bruzual & Chariot 2003). In Appendix A, we show that 
the commonly assumed forms of SFHs tend to system- 
atically bias the most probable estimates of the physical 
parameters. We make no assumptions here. Similarly, a 
common way to generate the SFH library is by Monte 
Carlo realizations of different SFHs. Instead, our library 
consists of a uniform, multi-dimensional grid of models 
without any prior assumption of the relevant physical 
parameters. This is necessary for exploring the complex 
SFHs of low-mass dIrr galaxies. When creating the SFH 
library, we allow the relative SFR among different age 
bins to logarithmically vary by an order of 1.5. 

We use the Chariot & Bruzual (in preparation) stellar 
population models, which have implemented several im- 
portant improvements compared to the Bruzual & Char- 
lot (2003) models. In particular, the new models include 



the recently improved treatment (Marigo et al. 2008) 
of the thermally pulsing asymptotic giant branch (TP- 
AGB), which at maximum (from a few hundred Myr to 
~ 1 Gyr for a single stellar population) could dominate 
the NIR emission. We adopt the stellar initial mass func- 
tion (IMF) parametrized by Chabrier (2003). We did not 
consider chemical evolution in the models, all the compo- 
nents of each model in the library have the same metal- 
licity. The metallicities of the population models are al- 
lowed to vary among 0.02, 0.2, and 0.4 Z0, which are 
adequate for our dIrr galaxies. 

To make use of the Ha data, the nebular hydrogen 
emission lines from H 11 regions are also obtained for each 
model in the library. CoUisionless Case B recombination 
for a 10,000 K gas is assumed. Using the recombination 
rates given by Hummer & Sorey (1987), the Ha luminos- 
ity is calculated from the UV stellar photons shortward 
of the Lyman limit. A systematic decline of the ratio of 
the integrated Ha-to-FUV flux in low luminosity dwarf 
galaxies has been shown (e.g. Sullivan et al. 2004; Meurer 
et al. 2009; Lee et al. 2009a). Three interpretations have 
been suggested in the literature, namely, the leakage of 
ionizing photons (e.g. Hunter et al. 2010), the stochas- 
ticity of massive SF at low SFR (e.g. Cervino & Luridi- 
ana 2004), and non- universal stellar IMF (especially the 
upper end; e.g. Pflamm-Altenburg, Weidner & Kroupa 
2009). These possibilities all imply that the observed Ha 
emission only gives a lower limit to the current SFR, 
under the assumption of Case B recombination and a 
well-populated universal IMF, which is the case for our 
modeling. Therefore, our Ha data are used to set a lower 
limit on the Ha fiux from the models. Considering the 
uncertainties in the observed Ha flux (HE04) , we set the 
lower limit to be 0.8xLHa,obs- 

We use the two-phase dust attenuation recipes devel- 
oped by Chariot & Fall (2000) as our standard dust ex- 
tinction law. The effective absorption is proportional to 
A^" '^. Chariot & Fall (2000) found a typical value of ~3 
for the ratio /i of the extinction experienced by the young 
stellar populations (< 10 Myr) to that by the older stellar 
populations. We adopt this value when creating our SFH 
library. The F-band dust extinction affecting the young 
stellar populations is allowed to vary linearly between 0.0 
and 0.5 magnitude. The gray extinction curve derived 
by Chariot & Fall (2000) is in agreement with recent re- 
sults for a large sample of nearby galaxies (Johnson et 
al. 2007). 

4.2. Determining the physical parameters 

Our final library consists of '--^ 4 x 10® different SFHs, 
which are created by allowing all the relevant parameters 
(i.e. dust extinction, metallicity, and relative SFR among 
different age bins) to vary uniformly (either logarithmi- 
cally or linearly) among physically reasonable ranges. For 
each SFH in the library, we integrate respectively back to 
the past 0.03 Gyr, 0.1 Gyr, 1 Gyr, and the galaxy lifetime 
(here we adopt a Hubble time of 13.7 Gyr) in order to get 
the SFR averaged over these different timescales. In the 
following we denote the SFR averaged over the past 0.03 
Gyr, 0.1 Gyr, 1 Gyr, and the whole lifetime as SFR0.03, 
SFRo.i, SFRi, and SFR13.7, respectively. The accumu- 
lated present-day stellar mass is also derived. From our 
experiments, the broadband SED modeling is not ex- 
pected to be very sensitive to the SFR averaged over the 
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past 0.03 Gyr (see the Appendix) and that averaged over 
the past 4 Gyr. So we concentrate on the interpretation 
of SFRo.i, SFRi, SFR13.7 and the stellar mass in this 
work. However, as mentioned above, our choice of six 
age bins guarantees that we have a relatively complete 
SFH library adequate for broadband SED modeling. 

We use the Bayesian technique (e.g. KaufFmann et al. 
2003; Kong et al. 2004) to find the most probable es- 
timates of parameters related to each of the observed 
broadband SEDs. Specifically, given an observed SED, 
we construct the probability density function (PDF) and 
the related cumulative distribution function (CDF) for 
each parameter defining our SFH library by calculating 
the likelihood exp(-x^/2) that an observed SED corre- 
sponds to each SFH model in our library. As described 
in the previous subsection, all the models with the emer- 
gent Ha luminosity (scaled by a factor determined from 
the weighted least squares) smaller than O.SxinQ.obs are 
excluded from the final determination of the parameters. 
The most probable value for a specific parameter is taken 
as the median of the corresponding PDF. The confidence 
interval is defined as the central 68% of the CDF. Simply 
adopting the single best-fitting value as the most prob- 
able estimate for the parameter is blind to degeneracies 
among relevant parameters. Our method is very robust, 
and the possible degeneracies are factored into the re- 
lated confidence intervals. In Appendix A, we thoroughly 
discuss the validity of these parameter estimates from the 
broadband SED modeling. 

4.3. Correcting the GALEX NUV data for Galactic 
extinction 

Before feeding the data into the SED modeling, we cor- 
rect the data for the foreground MW extinction (Schlegel, 
Finkbeiner & Davis 1998), adopting the Cardelli et al. 
(1989) extinction curve with Ry = Av/E{B -V) = 3.1. 
While there are nearly fixed, linear relations between 
Ax and E{B — V) for most bands, this is not the case 
for GALEX NUV (e.g. Wyder et al. 2007), due to the 
finite bandwidth (~ 1000 A) and the presence of the 
2175 A bump in the MW extinction curve. The ra- 
tio Afpj\i/E{B — V) is sensitive to both the emergent 
extra-galactic SED and the foreground MW reddening 
E{B — V). For example, the SED emergent from a con- 
stant SF and the SED from a 7 Gyr old single stellar 
population have A^xsm/ E{B — V)^ 8.2 and ^ 7.5 respec- 
tively for 0.05 mag of MW reddening. Also, for the SED 
emergent from a constant SF, Afpj\i/E{B — V) would be 
7.9 if the Galactic reddening was 0.5 mag. Therefore, 
unlike the other bands, the NUV data are not corrected 
for MW extinction with a fixed ratio of A-sxsm/EIB — V). 
Before comparing the data with each model in the li- 
brary, the NUV is corrected for the Galactic extinction 
using the ratio of A^xs\i/E{B — V) predetermined from the 
emergent SED of the model and the Galactic E{B — V) 
towards the specific galaxy. In other words, the MW ex- 
tinction correction for NUV is consistently applied before 
comparing each model with a given galaxy. 

5. RESULTS OF SED MODELING 

The SFH variations as a function of radius within the 
galaxies, determined from our SED modeling, are shown 
for each galaxy in the middle column of Figure [TJ Specif- 
ically, the SFR averaged over the past 0.03 Gyr, 0.1 Gyr, 



1 Gyr and the Hubble time are shown, respectively, as 
blue circles, green squares, brown triangles and red dia- 
monds. Stellar mass surface density profiles are shown 
in the third column of Figure [T] Table 3 gives the glob- 
ally integrated stellar mass, lifetime averaged SFR13.7, 
and ratios of SFRo.i to SFR13.7 and SFRi to SFR13.7 
for each galaxy. All the photometry with uncertainties 
a > 0.3 mag are excluded from the SED modeling. The 
global results are derived by integrating the relevant ra- 
dial profiles across the whole disk. 

5.1. Comparison with the stellar CMD-based SFHs 

Before proceeding to the interpretation of our SED 
modeling results, we compare our results with those de- 
rived from existing stellar CMD-based analysis. Half of 
our galaxies have stellar CMD-based SFHs obtained from 
HST data and available in the literature. Among these 
galaxies, 14 galaxies have CMD analysis for a signifi- 
cant part (the observing field covers more than 20% of 
the area enclosed by the Holmberg radius) of the stel- 
lar disk. The availability of GALEX FUV and Spitzer 
3.6 /im data makes our results for SFRo.i and M^ par- 
ticularly reliable. Therefore, we compare our results of 
the fraction of the SF accumulated during the recent 1 
Gyr fiGyr ((SFRi/SFRi3.7)/13.7) with that derived from 
CMD modeling. Among the 14 galaxies, the ficyr for 
11 galaxies (denoted as black circles in Figure |4]) are 
from Weisz et al. (2011), and those of the other 3 galax- 
ies (denoted as triangles in Figure [3|) are from Orban 
et al. (2008). From the comparison of these 14 galax- 
ies (Figure ^ , our globally integrated results from the 
SED-modeling are broadly consistent with those from the 
CMD-based analysis. 

Generally speaking, studying the resolved stellar popu- 
lations is the best way to derive a detailed SFH, provided 
the data are deep enough. We note that most of the above 
CMD analyses are carried out on relatively shallow data 
sets {Mi < 0). This means that the SFH at older times (> 

2 Gyr) is mainly constrained by the age-insensitive red 
giant branch (RGB). The RGB alone at best could pro- 
vide a good estimate of the total SF, from its formation 
to a few Gyr ago (Greggio 2002) to within a factor of ~ 
2. Nevertheless, the results from the current CMD-fitting 
are also sensitive to the modeling techniques used by dif- 
ferent works. For example, using the same HST data sets 
of WLM, Weisz et al. (2008) derived a two times higher 
fiGyr than Dolphin et al. (2005). Even more, adopting 
different stellar evolution libraries can lead to differences 
of at least a factor of 2 in the derived SFR (Gallart et 
al. 2005). Our library of SFHs are created with stellar 
evolution tracks combining those of Padova 1994 (see 
Bruzual & Chariot 2003) and Marigo et al. (2008) for 
the evolution of AGB, whereas the above CMD-modeling 
adopted the tracks of Padova 2000 (Girardi et al. 2002) 
and Marigo et al. (2008). This undoubtedly leads to some 
differences in the SFHs. We also noticed that the CMD- 
modeling works listed above adopted different IMFs from 
ours. Finally, the CMD-modeling often covers a smaller 
spatial area on the galaxy than our SED fitting. Given 
these uncertainties and differences, the results from our 
SED-modeling are in reasonable agreement with those 
from CMD-modeling. 

SFHs derived from broadband SED-modeling have 
been thought in the past to be biased towards younger. 
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luminous populations. The above comparison shows that 
the average SFR during the past ^ Gyr is well con- 
strained by modeling the broadband SED with a rel- 
atively complete library of SFHs. We point out that 
the availability of U and 5-band data is essential in 
SED-modeling because the two bands straddle the age- 
sensitive 4000 A break. 

5.2. Globally Integrated SFHs 

By integrating the relevant surface density profiles 
times the circumference over radius, we derived (in 
Table 3) the asymptotic total stellar masses M^, (See 

§ 5.4), SFR13.7, bo.lGyr (SFR0.1/SFR13.7), and biCyr 

(SFR1/SFR13.7). The galaxies are listed in order of in- 
creasing total baryonic mass in Table 3. A factor of 2-3 
(e.g. Hunter & Gallagher 1986; Kennicutt et al. 2005; 
Lee et al. 2009b) enhancement of current SF compared 
the averaged past rate has been used as one way to de- 
fine a starburst in the literature. Compared to SFR13.7, 
38% (13) of the sample galaxies have SFRq.i enhanced 
at least by a factor of 2, and 47% (16) have SFRg.i 
consistent with SFR13.7 within a factor of 2. For SFRi, 
12% (4 galaxies) have been enhanced at least by a fac- 
tor of 2, and 71% (24 galaxies) have SFRi consistent 
with SFR13 7 within a factor of 2. Among the 5 galaxies 
with SFR0.1/SFR13.7 < 0.5, 4 (LGS 3, DDO 216, DDO 
101, WLM) have baryonic mass (here approximated to be 
the stellar mass plus atomic gas mass, see below) smaller 
than 10^ Mq. Similarly, all 6 galaxies with SFR1/SFR13.7 
< 0.5 have baryonic mass smaller than 10* Mq. The total 
atomic gas mass was derived by multiplying the Hi gas 
mass by 1.34 to account for He. The H i masses used here 
were all from single dish H i emission line observations re- 
ported in the literature (see HE04 for the references), and 
the total stellar mass is derived from our SED modeling 
(see 

We plot ho iGyr and hicyr for the galaxies (Table 3) in 
Figure [5l The galaxies with baryonic mass smaller and 
higher than 10* M© are denoted as different symbols. A 
mass of 10* M© is about the median baryonic mass of our 
whole sample. There is almost no relationship between 
SFHs and total baryonic mass for our sample, except that 
all galaxies that exhibit extremely declining SFHs are 
relatively low mass systems. Among the whole sample, 
five galaxies (LGS 3, DDO 216, NGC 4163, DDO 101 
and WLM) have experienced a significant (>30%) decline 
in SF activity (both SFRq.i and SFRi) over the past ^ 
Gyr. We note that all of these five galaxies have baryonic 
masses smaller than 10* Mq, and, except for WLM, the 
(atomic) gas to stellar mass ratios are smaller than 1. 

The tidal indices (Table 1) of LGS 3 and DDO 216 
are larger than 1, which indicates that their evolution 
has likely been significantly influenced by environment. 
Most of the gas in these two galaxies could have been 
removed by either ram pressure stripping or tidal dis- 
turbance. In particular, based on the morphology of the 
Hi distribution, McConnachie et al. (2007) concluded 
that DDO 216 is currently undergoing ram pressure strip- 
ping caused by the intragroup/intracluster medium. The 
tidal index (Karachentsev et al. 2004) of a galaxy is de- 
fined as the maximum density enhancement caused by 
all neighboring galaxies, and higher values correspond to 
the likelihood of more significant tidal interaction with 



neighboring galaxies. This is consistent with the expec- 
tation that, lower mass systems (below a few times 10* 
M0) are more susceptible to significant influence of both 
the environment (Gunn & Gott 1972) and stellar feed- 
back (e.g. winds and supernova explosions (SNe)), which 
could even induce blow-out of interstellar gas from the 
galaxy (e.g. Mac Low & Ferrara 1999). 

A few galaxies deserve special mention. M81dwA, 
DDO 75, Haro 29 and NGC 1569 hold the largest ratio 
(> 5) of SFR0.1/SFR13.7 among our sample. The recent 
starburst in NGC 1569, which has the highest current 
SF intensity in our sample, may be ascribed to an in- 
teraction of some kind (e.g. Stil & Israel 1998; Muhle 
et al. 2005; Johnson et al. 2011). Similarly, the recent 
rise of the SFR in M81dwA may be due to the recent 
interaction within the M81 group, which has several dra- 
matic galaxy-galaxy interactions in progress. Haro 29 is 
a typical BCD galaxy undergoing a intense burst of SF 
in the central regions (Thuan & Martin 1981). DDO 75 
and Haro 29 both have negative tidal indices (-0.6 for 
DDO 75, -1 for Haro 29), indicating that the recent en- 
hancement of SF may be ascribed to internal processes. 
For DDO 154, Kennicutt & SkiUman (2001) reported a 
factor of ^ 2-4 times lower recent SF compared to the 
past based on Ha imaging, which is inconsistent with our 
results for SFRq.i. We interpret the disagreement as be- 
ing due to the fact that Ha is not a robust SFR indicator 
in LSB dwarf galaxies (e.g. Meurer et al. 2009; Lee et al. 
2009a; Hunter et al. 2010). 

5.3. Azimuthally-averaged Radial SFHs 

The sample as a whole exhibits a diversity of radial 
variations of SFHs (Figure [T]). However, some systematic 
trends do exist. 

Figures H] and [7] show the radial variations of 
log(I]sFRo.i/S*) and of log(I]sFRi/S*). The radius has been 
normalized by the V^-band disk scale length. The galaxies 
in Figures H] and [7] are plotted in order of increasing to- 
tal baryonic mass from the upper left panel to the lower 
right panel, and the galaxy names are listed in order 
of increasing baryonic mass in each panel. The dashed 
line in each panel marks a constant SFH over a Hubble 
time. Deep stellar CMD analysis of nearby dirr galax- 
ies (e.g. IC 1613: Skfllman et al. 2003; Leo A: Cole et 
al. 2007) suggests that most of the SF in dIrr galaxies 
may take place at intermediate ages, rather than persist- 
ing constantly over a Hubble time. In constrast, LGS 3 
(Hidalgo et al. 2011) consists of mostly old populations. 
Therefore, the data points lying above the dashed line 
in Figures [B] and [7] do not necessarily mean elevated SFR 
compared to the past. 

To quantify the radial trends of the SFH variation, we 
did linear least-squares fltting to the radial variations 
of log(i;sFRo.i/S*) and log(i;sFRi/S*). If there is an ob- 
vious break in the stellar mass surface density profiles 
(Figure [T] and Table 3), we fit the slope for the inner 
and the outer parts separately. The GALEX data (thus 
the measure of recent SFR) of several galaxies (LGS 3, 
DDO 210, DDO 216 and NGC 1569) have high S/N ra- 
tios only in the inner disk, in this case we only fit the 
inner disk. We also fit the relevant ratios as a function 
of normalized (by the F-band scale length) radius. The 
slopes (or gradients) from the fitting are listed in Ta- 
ble 4. The weighted averages of the relevant slopes for 
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the galaxies below and above a baryonic mass 10^ Mq 
are also listed in Table 4. The slopes are further plot- 
ted in Figures [8] and ID Also shown in Figures |8] and |9] 
are the weighted averages {large diamonds) of the whole 
sample. Negative slope indicates that the SFR averaged 
over the relevant timescale becomes more centrally con- 
centrated compared to the past. As is shown, on average, 
our sample galaxies have more negative radial variations 
of log(EsFR„ j/E*) and log(EsFRi/E^) in the outer disks 
than the inner disks, and log(I]sFRo i/E*) has steeper ra- 
dial declining than log(EsFRj/E*). In the inner disks, 44% 
(15 galaxies) of the sample show negative radial slopes 
for both log(SsFR„,i/E^) and log(EsFRi/Ei,), and 9 of the 
15 galaxies have M^ary < 10® Mq. A few galaxies with 
Mbary > 10® Mq (i.e. DDO 165, DDO 63, DDO 154, 
DDO 168, NGC 2366 and NGC 4214) have almost flat 
slopes (|slope| < 0.2) in the inner disks. In the outer 
disks, 80% (27 galaxies) have negative slopes for both 
log(i;sFRo.i/E*) and log(EsFRi/E^), and 14 of the 27 galax- 
ies have Mbary < 10® Mq. Particularly, 35% (12 galaxies, 
7 with Mbary ^ 10® Mq) of the galaxies exhibit nega- 
tive slopes across the observed disk. In NGC 4163 and 
NGC 3738, the recent SF has been suppressed by more 
than an order of magnitude in the outer regions, even 
though these two galaxies have almost constant ratios of 
log(i;sFRo.i/E^) and log(SsFRi/E*) at large radii. These 
results are in line with Figure [2l which shows the com- 
parison between the scale length of FUV, B and 3.6 /im. 

Among the whole sample, the weighted averages of the 
slopes of A(log(EsFR„.i/E*))/A(R) are 0.10±0.42 and - 
0.64±0.48 for the inner and the outer disks, respectively. 
The weighted averages of A(log(I]sFRi/Ei))/A(R) are - 
0.05±0.37 and -0.28±0.39 for the inner and the outer 
disks, respectively. Furthermore, as listed in Table 4, 
the inner disks tend to have much shallower radial slopes 
of A(log(EsFR/E^))/A(R) than the outer disks. The SF 
activity in dirr galaxies was suggested to be a random 
percolating process across the stellar disk (e.g. van Zee 
2001). According to our results, this randomly perco- 
lating scenario may be only appropriate for the inner 
disks of some galaxies (especially those relatively mas- 
sive galaxies). As far as the whole disk of the low mass 
dIrr galaxy is concerned, the star-forming disk has been 
shrinking from the outer part. 

5.4. Stellar Mass Surface Density Profiles 

The stellar mass surface density E^ profiles of the 
galaxies are shown in the last column of Figure 1. The 
stellar mass surface density profiles are well fitted with a 
piece-wise exponential function, which is overplotted on 
the profiles. For some relatively faint galaxies, the high 
quality IRAC 3.6 fim data reaches further radii than the 
other bands. In this case, considering the relative insen- 
sitivity of the Mj,/Lnir to the underlying stellar popula- 
tions (Bell & de Jong 2001), we extend stellar mass sur- 
face density profiles derived from our multi-band fitting 
to the radius reached by the 3.6 /im surface photometry, 
adopting the M^,/ L^,e^rn ratio from the radius where our 
multi-band SED fitting ends. The 3.6 /im photometry for 
some galaxies with large angular size does not reach as 
far as the other bands due to the limited field of view 
of the Spitzer observations. In this case, we constrained 
the stellar mass profiles at the outer part by modeling the 



UBV{&nd G^LSXFUV/NUV if available) data. We also 
derived the asymptotic total stellar mass by applying the 
"growth curve" method to the accumulated stellar mass 
surface density profiles (e.g. Mufioz-Mateos et al. 2009). 
Briefly, we obtained the radial gradient of the accumu- 
lated stellar mass surface density at each observed ra- 
dius, then the accumulated stellar mass surface density 
was linearly fit as a function of the gradient and the y- 
intercept (zero gradient) of the fit was adopted as the 
asymptotic stellar mass. Table 3 lists the central stel- 
lar mass density E^-^center extrapolated from the inner 
exponential disk, stellar mass concentration index C31, 
the inner/outer scale lengths of the stellar mass density 
profiles, and the asymptotic total stellar mass for each 
galaxy. C31 is defined here as the ratio of the radii that 
encompass 75% and 25% of the total stellar mass (de 
Vaucouleurs 1977). 

Like the radial light profiles, the stellar mass surface 
density profiles of dIrr galaxies usually can be described 
by either a single or piece-wise exponential profile over 
most of the disk. The most striking feature about our 
galaxies is that most of them have broken stellar mass 
surface density profiles. Specifically, ~ 80% (27) of the 
galaxies exhibit obvious down-bending profiles, in the 
sense that the outer disks have a steeper mass profile than 
the inner disks. Whereas three (DDO 70, Haro 29, NGC 
1569) galaxies have obvious up-bending (flatter outer 
disk) profiles, in the sense that the outer disks exhibit 
shallower profiles than the inner disks. Both Haro 29 and 
NGC 1569 have experienced intense, extended starbursts 
in their inner regions. It has been shown that down- 
bending surface brightness profiles are very common in 
spiral galaxies (e.g. Ferguson & Clarke 2001; Pohlen & 
Trujillo 2006). Nevertheless, Bakos, Trujillo & Pohlen 
(2008) found a pure exponential stellar mass profile for 
the spiral galaxies with down-bending surface brightness 
profiles. This is different from our dIrr galaxies which 
tend to have both down-bending stellar mass profiles and 
down-bending surface brightness profiles. Amorin et al. 
(2007) found that the LSB steUar hosts in BCD galax- 
ies have near-exponential profiles. From our results for 
BCD galaxies, not only the LSB stellar hosts, but also 
the inner starburst regions have exponential stellar mass 
surface density profiles, although the scale length is dif- 
ferent for the inner starburst region and the LSB stellar 
host (Figure [T]). 

The present-day stellar mass surface density profile of 
dIrr galaxies should be primarily determined (see § 6.1.1) 
by the radial variations of SFH. As is shown above (Ta- 
ble 4), on average, the SF in the outer disks has been 
decreasing more significantly both over time (Figures [6] 
and[7]) and over radius (Figures[8]and[9]) compared to the 
inner disks, which would naturally lead to down-bending 
stellar mass surface density profile, seen in most of our 
galaxies. This scenario of outside-in depression of SF is 
different from what may be happening in typical BCD 
(e.g. Haro 29, NGC 1569, see FigureH]) galaxies. The cen- 
tral starburst of BCD galaxies results in faster buildup 
of the inner stellar disk than the outer part, which could 
lead to up-bending stellar mass surface density profiles. 

As can be seen, a local enhancement of the recent SF 
may affect the stellar mass density profile in these low 
mass dIrr galaxies. For example, in NGC 2366, the su- 
pergiant Hii complex (NGC 2363, Youngblood & Hunter 
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1999) ~ 1 kpc away from the center corresponds to the 
remarkable spike in the stellar mass surface density pro- 
file. The concentration index of the stellar mass surface 
density profile is generally close to, but lower than, that 
of a single exponential (C31 = 2.81). 

Lower mass systems tend to have lower stellar mass 
concentrations (Tabic 3). However, the BCD galaxy Haro 
29 is as centrally concentrated (C31 = 5.3) as a de Vau- 
couleurs r^/^ bulge profile, which suggests that the cen- 
tralized starburst has lasted an extended period of time. 
Assuming the SFRq.i in the central regions is represen- 
tative of the recent episode of central starburst, and the 
central stellar mass is dominantly contributed by the 
burst, the starburst would have lasted more than 0.6 Gyr. 

6. DISCUSSIONS 

6.1. The disk assembly mode of dirr galaxies: outside-in 

The "inside-out" growth mode has long been suggested 
for the formation of galaxy disks (e.g. Larson 1976; Chi- 
appini et al. 1997; Mo, Mao & White 1998; Naab & Os- 
triker 2006). The "inside-out" scenario reproduces many 
observations of spiral galaxies, such as the radial gra- 
dients of both broad-band colors (bluer outward) and 
metallicity (lower outward) (dc Jong 1996; Bell & de 
Jong 2000; MacArthur et al. 2004; Wang et al. 2011), and 
the extended UV emission discovered in outer spiral disks 
(e.g. Thilker ct al. 2007; Boissier et al. 2008). Muhoz- 
Mateos et al. (2007) found a moderate inside-out disk 
formation by studying the radial profiles of G ALEX FW 
and 2MASS Ks for a sample of relatively face-on nearby 
spiral galaxies. Recently, the "inside-out" growth mode 
has been confirmed for M 33 (Williams et al. 2009) and 
NGC 300 (Gogarten et al. 2010) from analysis of CMDs 
obtained with HST. For M 33, it was shown (Williams 
ct al. 2009; Barker ct al. 2011) that the "inside-out" sce- 
nario only applies to the inner disk, and the region be- 
yond the surface brightness break exhibits positive age 
gradient. 

However, unlike luminous spiral galaxies, late- type dIrr 
galaxies exhibit a variety of behaviors in terms of color 
profiles (Kormendy & Djorgovski 1989; TuUy et al. 1996; 
Jansen et al. 2000; Taylor et al. 2005; HE06). Neverthe- 
less, Tully et al. (1996) found that galaxies fainter than 
Mb ~ -17 become redder with radius in the Ursa Major 
cluster. Later, with a larger sample, Jansen et al. (2000) 
found that the reddening (with radius) trend may occur 
at a fainter absolute magnitude than that found by Tully 
ct al. (1996). Of a sample of 94 dIrr galaxies studied by 
HE06, 64% of the galaxies with a gradient in B — V be- 
come redder with radius. In particular, all the galaxies 
with Mb > -14 become redder at least in one color (ei- 
ther B — V or U — B) at outer radii. Recently, Tortora 
et al. (2010) analyzed the optical color gradients of ~ 
50,000 nearby Sloan Digital Sky Survey galaxies. They 
found a good correspondence between the color gradients 
and the stellar mass of the galaxies. In particular, below 
a stellar mass M^, < 10^ '' Mq, the color gradient slopes 
become positive, which were mainly attributed to the 
positive metallicity gradients. The lowest mass galaxies 
in the sample of Tortora et al. (2010) have - 10^ ^ 
Mq, compared to the median stellar mass ^ 10^'^ Mq 
of our sample. All the above studies are based on optical 
broadband data alone, which makes the interpretation of 



the radial color gradients ambiguous, d\ie to the degen- 
eracy between age, metallicity, and extinction. 

Additionally, from stellar CMD analysis of two dif- 
ferent (HST/WFPC2) fields in IC1613, Skillman et al. 
(2003) (also Bernard et al. 2007) found that the SF ac- 
tivity in the outer field has been significantly depressed 
during the last Gyr, which is in contrast to the "inside- 
out" growth scenario. Similarly, by analyzing resolved 
stellar populations, the characteristics of more spatially 
extended older populations have been reported both 
for transition-type dwarf galaxies (e.g. DDO 210: Mc- 
Connachic et al. 2006; Phoenix: Hidalgo et al. 2009) and 
dwarf spheroidal (dSph) galaxies (e.g. Sculptor dSph: 
Tolstoy ct al. 2004; Fornax dSph: Battaglia ct al. 2006). 
A clear trend of "outside-in" quenching of recent SF was 
also found (Gallart et al. 2008; Indu & Subramaniam 
2011) in the Large Magellanic Cloud, which has a mass 
comparable to that of the most massive galaxy (i.e. NGC 
1156) in our sample. 

These studies are in agreement with our findings for 
our sample of dIrr galaxies. We point out that our anal- 
ysis here is only sensitive to the SF during the recent ~ 
Gyr and the SF averaged over the whole lifetime. The 
radial variations of SF at the intermediate ages could 
be different from the recent past. For instance, in the 
late-type spiral galaxy NGC 2976, Williams et al. (2010) 
found a deficit of stellar populations younger than ~ Gyr 
beyond the break of the disk surface brightness profile, 
with similar ancient stellar populations at all radii. Nev- 
ertheless, since the suppression of SF in the outer disk 
was found in the majority of our sample galaxies, the 
trend we found here must reflect the evolutionary process 
of the stellar disks of low-mass dIrr galaxies in general. 
In the following, we discuss possible interpretations of 
the observed radial variations of the stellar populations 
in dIrr galaxies. 

6.1.1. In situ SF vs. Secular redistribution? 

The present-day radial distribution of the stellar pop- 
ulations could be ascribed either to in situ SF or to 
a stellar redistribution process. For example, disk non- 
axisymmctric instabilities (e.g. spiral arms, bars) arc ca- 
pable of a large-scale redistribution of stars. In spiral 
galaxies, resonant scattering with transient spiral arms 
can lead to substantial radial steUar migration (e.g. Sell- 
wood & Binney 2002; Roskar et al. 2008), which compli- 
cates determination of the in situ SFH. Likewise, the 
existence of bars can, besides inducing strong gas in- 
flows, drive substantial redistribution of the stellar disk, 
civen increasing the disk scale length (e.g. Hohl 1971; De- 
battista et al. 2006). However, dIrr galaxies usually lack 
large-scale instabilities (H98), which is a requisite of gen- 
crating spiral arms and (probably) bar instabilities (e.g. 
Lin & Pringlc 1987; Mihos, McGaugh, & de Blok 1997). 

Nevertheless, it has been shown that a slowly rising 
rotation curve naturally leads to the alignment of elon- 
gated orbits, and thus is more prone to bar formation 
(Lynden-Bell 1979). Even more, once the bar is made, it 
may exist for many dynamical times in a situation of solid 
body rotation, which is the case over a large part of the 
disk in most dIrr galaxies. Nevertheless, the bar growth 
of late-type galaxies may stop at a very early stage be- 
cause the co-rotation radius moves out of the main disk 
(Combes & Elmegreen 1993). There exists a strong cor- 
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relation between the strength of the bar and the central 
density (e.g. Elmegreen et al. 2007; Sheth et al. 2008). 
dirr galaxies usually lack a central stellar mass excess, 
which implies that bar instabilities, if they exist, are very 
inefficient in transferring angular momentum across the 
disk of dIrr galaxies. Of our sample, 11 galaxies (32%) 
show evidence for a bar based on ellipse fitting of opti- 
cal isophotes (HE06). Figures [8] and [9] suggest that there 
is no preference for barred galaxies exhibiting larger or 
smaller radial gradients of log(S]sFR/S*). We note that 
identifying bars in dIrr galaxies is not as straightforward 
as in luminous spiral galaxies. 

6.1.2. External influences 

External factors, such as interactions (both with the 
intragroup/intracluster medium and with neighboring 
galaxies) and the cosmic UV background (e.g. Gnedin 
2000), can also lead to more centrally concentrated re- 
cent SF or depressed SF in outer disks. The tidal in- 
dex (Table 1) suggests that less than 5 of the galaxies 
studied here could be noticeably affected (© ^ 1) by 
neighboring galaxies. The slopes of the radial variations 
of log(i;sFRo.i/S*) and log(EsFRi/S^) for the outer disks 
are plotted against the tidal index in Figure [10] (There 
is no calculation of 9 for DDO 101 in Karachentsev et 
al. (2004), so we arbitrarily set 8 as when creating 
the plot). There is no obvious correlation between the 
radial variations of the SFHs and the environment. How- 
ever, most galaxies fall between an upper envelope and 
a lower envelope, as indicated by the dashed lines in the 
plot. NGC 6822 has such a large angular size projected 
on the sky that the observations only covered the cen- 
tral regions. As mentioned above, the SF of dIrr galax- 
ies usually peaks in the circumnuclear regions, therefore 
the slopes measured here for NGC 6822 do not reflect 
the SFH variations across the main star-forming disk. In 
fact, based on HST images of five fields in NGC 6822, 
Wyder (2001) found a higher ratio of the recent SFR to 
the average past rate in the inner bar regions than the 
outer regions, consistent with the trend we found in this 
study. The galaxies with larger Q prefer lower (zero or 
negative) slopes of the radial variations of log(I]sFRo J^*) 
and log(i;sFRi/S*). 

Two of the most obvious environmental effects stud- 
ied in the literature are ram pressure stripping of 
the gas component caused by a hot gaseous intra- 
group/intracluster medium (Gunn & Gott 1972; Lin 
& Faber 1983) and tidal disturbance from neighboring 
galaxies (e.g. Mayer et al. 2001). By studying the dy- 
namical properties of the dwarf galaxies in the Fornax 
cluster, Drinkwater et al. (2001) show that the dwarf 
galaxies are still falling into the cluster, and the frac- 
tion of dwarfs with active SF drops rapidly towards the 
cluster center. This is a clear indication of the environ- 
mental effect on the evolution of dwarf galaxies. Simi- 
larly, Pustilnik et al. (2002) found that the blue compact 
galaxies in higher density environment have on average 
less H I which could be attributed to either ram pressure 
stripping of the gas or tidal interactions. 

Presumably, both ram pressure stripping and tidal dis- 
turbance should be more effective on lower mass galax- 
ies. The outer disk, especially in relatively low mass sys- 
tems where the gravitational well and the gas surface 
density are relatively low, is more prone to gas removal 



due to ram pressure. Furthermore, significant tidal forces 
from neighboring galaxies could induce a bar instability 
and ensuing gas inflows towards the inner regions (e.g. 
Mayer et al. 2001), which leads to a more centrally con- 
centrated SF. In our sample, lower mass galaxies tend 
to have more centrally concentrated SF than relatively 
higher mass galaxies, which is in general agreement with 
the expectation of the above two environmental effects. 

By studying the H I distribution of spiral galaxies in 
the Virgo cluster, Cayatte et al. (1994) found that some 
galaxies close to the cluster core have normal gas densi- 
ties (compared to field counterparts) in the inner disks 
but show a strong decline or cutoff of the Hi intensity 
starting within the half light radius of the stellar disk. 
These galaxies should have been undergoing remarkable 
ram pressure stripping. Actually, the radial H i distribu- 
tion (Hunter et al. 2011, in preparation) for some of our 
galaxies do show a strong decline starting well within the 
stellar disk. Numerical simulations (e.g. Vollmer 2003; 
Kapferer et al. 2008) suggest that both ram pressure 
stripping and tidal disturbance are needed to reproduce 
the observations. FigurefTUlsuggests that the tidal distur- 
bance from neighboring galaxies should play some role in 
the inward movement of SF activity. 

The presence of a UV background can heat the gas 
in low mass halos, which prevents SF at least in the 
outer disk where the gas is not dense enough for self- 
shielding (Susa & Umemura 2004). By solving the radia- 
tive transfer equation for the diffuse UV background in 
a pre-galactic cloud, Tajiri & Umemura (1998) showed 
that, above a critical number density ricrit ^ 1-4 x 10^^ 
cm-3 (Nhi ~ 1.3 X 10^^ for a gas disk - 0.3 kpc thick), 
which is almost independent of the UV background in- 
tensity, self-shielding of the ISM against the UV back- 
ground is prominent. Based on our H i emission line maps 
of these galaxies, and according to this critical density, 
the typical gas column density (~ 10^° cm~^) across the 
observed stellar disks of our galaxies should be sufficient 
for them to be self-shielded. Therefore, the cosmic UV 
background is unlikely to be the main driver of the cen- 
trally concentrated SF. 

6.1.3. Regulation of the SFR through stellar-feedback 

Depressed SF in outer dwarf disks could be attributed 
to self-regulation due to stellar feedback. Feedback from 
both SNe and stellar winds has been shown to be a ma- 
jor factor shaping the evolution of dwarf galaxies due 
to their shallow potential wells (e.g. Dekel & Silk 1986; 
Mori, Ferrara & Madau 2002; Mashchenko, Wadsley, & 
Couchman 2008). From the numerical simulations, the 
SNe/stellar winds-regulated evolution of isolated dwarf 
galaxies tends to be characterized by an episodic SFH 
(e.g. Chiosi & Carraro 2002; Stinson et al. 2007). Such 
SFHs have been suggested to explain the anomalously 
blue colors and low metallicities of dIrr galaxies (e.g. 
Searle, Sargent, & Bagnuolo 1973). 

In this SNe/stellar winds-regulated evolution scenario, 
the SN explosions and stellar winds from massive stars 
following one episode of SF trigger gas heating and ex- 
pansion, and some gas may be blown away from the 
galaxy. Subsequent cooling allows the remaining gas to 
sink deeper into the center and start another cycle of 
SF. The outside-in formation scenario has been suggested 
for low mass elliptical galaxies (Martinelli, Matteucci & 
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Colafrancesco 1998; Pipino & Matteucci 2004) to ex- 
plain the positive gradients of [Mg/Fe] and the blue col- 
ors in the center of dwarf elhptical galaxies. The outer 
regions, where the gravitational potential well is rela- 
tively low, are more susceptible to the galactic wind (thus 
blow-out/blow-away) compared to the inner regions. The 
blow-out/blow-away of the high angular momentum gas 
(compared to the inner regions) in the outer part of the 
dwarf disk results in a remarkable net loss of angular 
momentum, and thus the shrinking of the gas disk. This 
leads to a depression of SF in the outskirts and more pro- 
longed SF in the central regions. Recent simulations (Val- 
cke et al. 2008) indeed found that, for a self-regulated, 
SNe-driven evolution of the low mass (Mbary ^ 10^ M©) 
dwarf galaxies, the SF becomes more and more centrally 
concentrated with time, which is in agreement with what 
we see here. We point out that this scenario of outside- 
in shrinking may not apply to the early universe when 
the galaxy just started assembling and the potential was 
shallow (susceptible to stellar feedback) even in the inner 
region. Governato et al. (2010) found that the outflows 
resulted from stellar feedback may selectively remove low 
angular momentum gas, especially at z > 1. 

Recent high-resolution hydrodynamical simulations in- 
cluding cooling, SN feedback and the UV background 
radiation (Sawala et al. 2010) suggest that, while the 
combined effects of SN feedback and UV heating are nec- 
essary to reproduce the observations, SN feedback is the 
key factor in determining the evolution of low mass dwarf 
galaxies. The cosmic UV radiation has almost no effect 
on galaxy evolution if the SN feedback was ignored in 
the simulations. 

6.1.4. Pressure support in the gas disk 

Low-mass dirr galaxies are dynamically hot systems 
compared to luminous spiral galaxies. The (turbulent) 
gas velocity dispersion within dwarf galaxies can be com- 
parable to the ordered rotational velocity, which means 
the gas pressure may provide significant radial support 
to the rotating gas disk (Stinson et al. 2009; Dalcanton 
& Stilp 2010). Since the gas density usually declines with 
radius, the pressure gradient, together with the centrifu- 
gal acceleration, will counteract the gravity. A higher ef- 
ficiency of SF in the inner regions then leads to quicker 
consumption of gas, which results in an initial reduction 
of the gas density there. After the enhanced turbulence 
from SF subsides, the pressure in the center will drop 
too, and more gas will flow in so that pressure gradients 
and centrifugal force balance gravity again. Thus some 
shrinking of the gas disk might be possible, within the 
limits of angular momentum conservation. In their simu- 
lations of isolated low mass dwarf galaxies, Stinson et al. 
(2009) found a positive age gradient with radius after gas 
contraction, stellar migration, and supernova-triggered 
SF, with the latter two effects being more important for 
lower mass galaxies. Dalcanton & Stilp (2010) further 
demonstrated that the impact of pressure support alone 
on the evolution of gas disk may be weak. 

6.2. The relationship between Ssfrj and SFH 

H98 found that the radial profiles of current SF activity 
correlate better with those of the older stellar emission 
(e.g. ]^-band) than with the atomic gas in dIrr galaxies. 



Since star-forming molecular clouds are forming from the 
atomic gas, this lack of a correlation between SF and gas 
is unexpected. Similarly, Ryder & Dopita (1994; see also 
James, Bretherton, & Knapen 2009) discovered a good 
correlation between /-band surface brightness (a proxy 
for the mass dominant old stellar populations) and Ha 
for a sample of nearby spiral galaxies. Recently, a good 
correspondence between the surface brightness in SDSS 
r-band and in FUV was also shown for a sample of LSB 
galaxies (Wyder et al. 2009). If the optical broadband 
emission traces the underlying mass-dominant old stellar 
populations, the above correlations would be equivalent 
to a relationship between SF and the stellar mass surface 
density. Is this correlation between Ssfr and causal 
or casual? 

In spiral galaxies, the stellar component usually dom- 
inates the inner disk in terms of both surface density 
and volume density, so we might expect that the stel- 
lar component plays an important role in regulating the 
ISM and thus the SF process (e.g. Wong & Blitz; Blitz & 
Rosolowsky 2006; Shi et al. 2011). However, that is not 
necessarily the case in dIrr galaxies, where the atomic gas 
usually dominates the baryonic component. In Figure [TTl 
SsFRo 1 is plotted against E^, for the whole sample. The 
data points are extracted from the azimuthal averages 
of the relevant surface density profiles. Although there is 
rough correspondence between Ssfro i and E^,, the scatter 
is substantial. The Spearman rank correlation coefficient 
between E* and Esfro i is ^ 0.6. 

As is shown in Figure |6l the inner regions of our 
galaxies tend to have higher ratios of Esfrq i to E^, than 
outer parts, which is especially significant for relatively 
lower mass systems. In other words, Esfro i usually has 
a steeper radial dependence than E^,. This is consistent 
with our finding that shorter wavelengths tend to have 
shorter scale lengths than longer wavelengths (Figure [2]). 
Within the inner disks, our sample galaxies, on average, 
have almost zero radial gradients of Esfr/E^, (Figures [S] 
and [9]). Therefore, the correlation between Esfr and E^ 
is only good within the inner disk. 

Recently, Leroy et al. (2008) found that, in dIrr galax- 
ies and outer disks of spiral galaxies, the ratio of SFR 
surface density to atomic gas surface density exhibits an 
almost linear relationship to the stellar mass surface den- 
sity. Given the narrow range of the H i surface density (5- 
10 Mq pc~^, Leroy et al. 2008) within the optical disk, 
the above correlation actually reflects the correlation be- 
tween the SFR and the stellar mass surface density. We 
note that all of the dIrr galaxies studied by Leroy et al. 
(2008) have baryonic masses well above 10* M©. For the 
four galaxies in common with the sample of Leroy et 
al. (i.e. DDO 154, DDO 63, DDO 50 and NGC 4214), 
except for DDO 50, all the other three have shallow 
(|A(log(EsFRo.i/E^))/A(R)| < 0.2, Table 4) radial slopes 
of log(EsFRo.i/Ei,) across the inner disks, which extend 
to more than ^ 3 disk scale lengths, in agreement with 
Leroy et al. (2008). Nevertheless, some galaxies (DDO 
52, NGC 1569, NGC 3738, and NGC 1156) in our sam- 
ple with baryonic masses larger than 10* M© show obvi- 
ous gradients (either positive or negative, see Table 4) of 
log(EsFRo.i/E^) across the whole disks. 

The SFE in dIrr galaxies generally is very low (e.g. van 
Zee 2001; Lee et al. 2011), as is evidenced by high gas- 
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to-stcUar mass ratios (e.g. sec Table 1). Therefore, the 
stellar disk may have always been a small perturbation 
on the gas-dominated baryonic disk. The constant aver- 
age ratios of Ssfro.i to and Ssfri to seen in the 
inner disks of dirr galaxies could imply that the SFE at 
a given radius in the inner disk may have been nearly 
the same for most of the galaxy's lifetime. This would 
naturally lead to similar radial profiles of the SFR to the 
accumulated stellar mass. In agreement with this evolu- 
tionary picture, Elmegreen & Hunter (2006) considered 
multi-component SF triggering processes, including tur- 
bulent compression, to interpret SF in outer spiral disks 
and in dIrr galaxies. They suggested that the SF pro- 
file should be about the same as the accumulated stellar 
mass profile when the baryonic component of the disk is 
still dominated by atomic gas. In this scenario, the stellar 
mass distribution naturally follows the inefficient, patchy 
SF. On the other hand, for the outer disks, the SF pro- 
cess has been significantly affected by stellar feedback 
(e.g. SNe) and the external disturbance. As discussed 
in the previous section, feedback-regulated evolution or 
environmental effects may result in shrinking of the star- 
forming gas disk toward the inner part, and thus more 
and more centrally concentrated SF, without a remark- 
able correlation between SFR and stellar mass. We no- 
ticed that "inner disk" of galaxies with Mbary >10^ M.q 
could extend over more than ~ 3 disk scale lengths. 

7. SUMMARY 

We present multi-band (GALEX FUV/NUV, UBV, 
Ha and Spitzer 3.6 fiia) surface brightness profiles for 
a representative sample of nearby dIrr galaxies. The me- 
dian of the stellar mass and of the total baryonic mass 
of our sample galaxies are ~ 10^'^ Mq and 10® Mq, re- 
spectively. By modeling the multi-band SEDs with a rel- 
atively complete library of SFHs, we constrained: the 
radial variations of the SFR averaged over the recent 0.1 
Gyr, SFRo.i, the SFR averaged over the recent 1 Gyr, 
SFRi, and the stellar mass surface density profiles. To 
summarize the results: 

1. We show that, with a relatively complete library of 
model SFHs, the averaged SFR during the past 1 
Gyr can be well constrained by modeling integrated 
multi-band data. The exponential-plus-burst SFH 
library commonly used for spirals tends to sys- 
tematically overestimate recent SF of dIrr galax- 
ies, which are characterized as having much more 
complex SFHs than the larger spiral galaxies. 

2. For ~ 80% (27) of the dIrr galaxies studied in this 

work, and all (15) with baryonic mass less than 10® 
M0, surface brightness profiles determined from 
shorter-wavelength passbands have shorter expo- 
nential disk scale lengths than those based on pass- 
bands at longer wavelengths at least in the outer 
disks. This trend suggests that the star-forming 
stellar disk may be shrinking for most dIrr galax- 
ies. This produces an "outside-in" disk formation 
scenario for dIrr galaxies. 

3. The recent SF activity of ^ 38% of our galaxies 
has been enhanced by a factor of at least 2 com- 
pared to the lifetime averaged SFR. Even though 
no obvious correlation between global SFHs and 



total baryonic mass was found in our sample, all 
of the galaxies with significantly declining global 
SF activity have baryonic masses smaller than 10® 
Mq. This may be ascribed to the susceptibility of 
low mass dwarf galaxies to significant influence by 
both internal and external processes. 

4. Consistent with the "outside-in" scenario, 80% of 
our sample galaxies have negative slopes of radial 
variations of log(i;sFRo.i/S*) and of log(EsFRi/S*) 
at least in the outer disks (beyond the breaks of 
stellar mass surface density profiles). In the in- 
ner disks, especially for the galaxies with Mbary 
>10® M0, the slopes are much shallower (consis- 
tent with zero slopes on average) than that in the 
outer disks. Both internal (stellar feedback, gas 
pressure support) and external (ram pressure strip- 
ping of gas and tidal disturbance) processes are dis- 
cussed as the possible explanation for the "outside- 
in" shrinking of the star-forming disk. 

5. The radial distribution of the stellar mass density 
of dIrr galaxies is well described as a single or, 
for the majority, a piece-wise exponential profile. 
In particular, the majority (^ 80%) of our sample 
galaxies exhibit steeper outer (down-bending) stel- 
lar mass profiles. Since the obiter disks on average 
have much steeper radial declining Ssfr/Sj, than 
the inner regions, we interpret the down-bending 
stellar mass surface density profiles as a natural 
result of a gradual shrinking of the star-forming 
disks. A few galaxies exhibit much steeper stellar 
mass profiles in the inner disk than the outer disk, 
which may be caused by an intense starburst in 
their inner regions. Spiral galaxies with double ex- 
ponential surface brightness profiles and inside-out 
star formation may require a different explanation 
for the steep outer profiles. In spirals, the outer 
disks are sometimes younger than the inner disks, 
while in dwarfs, the outer disks are usually older 
than the inner disks. Perhaps spirals are still ac- 
creting gas in their outer parts, whereas the dwarfs 
have lost the gas in their outer parts and have other 
dynamical effects not present in larger galaxies. 

6. The previously found correlation between surface 
density of the stellar mass and recent SFR is only 
remarkable in the inner disks, which could extend 
over more than 3 disk scale lengths for the rela- 
tively massive {Mbary ^ 10® Mq) dIrr galaxies. The 
correlation can be explained by inefficient, constant 
SF at a given radius during most of the galaxy's 
lifetime. For the outer disks of most (80% for the 
present sample) dIrr galaxies, on the other hand, 
the correlation between SFR and stellar mass sur- 
face densities is not strong, due to effects of inter- 
nal feedback, turbulent gas pressure support, and 
external disturbances. 
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APPENDIX 
DETAILS CONCERNING SED MODELING 
Recovery of physical parameters 

Here we explore the rehabihty of multi-band SED modeling in recovering the rough SFH of dirr galaxies. In these 
tests, we first logarithmically divided the Hubble time (13.7 Gyr) into 14 independent age bins, and took the SF as 
constant within each age bin. Then a series of 3000 mock SFHs were generated by randomly changing the SFR for 
each independent age bin. The internal extinction was fixed at Ay = 0.1 when generating the SEDs related to each 
mock SFH. We then modeled the SEDs (GALEX FUV/NUV, U, B, V, 3.6 /^m) related to the mock SFHs with our 
library of 6-component superposition (6-C) SFHs used in this work. The calibration uncertainties of the real data were 
taken as the photometric uncertainties of the SED related to each mock SFH. 

Figure [12] presents a comparison of the modeling results, including M^,, SFRi, SFRq.i and SFRg.oa, and relevant 
integrated quantities of the mock input SFHs. We fit each histogram in Figure [12] with a gaussian curve, which is 
overplotted as a red solid line in each panel. The standard deviation [a) and mean value (/i) of the gaussian curve is 
also indicated. As can be seen, except for SFR0.03, all the parameters are well recovered in our SED modeling, with 
acceptable uncertainties related to the varying SFHs. Figure [13] further shows ratios of SFR averaged over different 
timescales. Again, except for SFR0.03, the modeling with the 6-C library recovers the SFHs with negligible, if any, 
bias. 

Uncertainty in converting the GALEX FUV luminosity to SFR 

FUV has been extensively used to estimate the recent (~ 0.1 Gyr) SFR in galaxies (e.g. Kennicutt 1998; Salim et 
al. 2007; Hunter et al. 2010). The most commonly used formula (e.g. Kennicutt 1998) assumes constant SF during the 
past 0.1 Gyr, which is comparable to a galactic dynamical timescale. The assumption of more or less constant recent 
SF is generally true for a high luminosity galaxy as a whole, as is evidenced by the consistency between the Ha-based 
SFR and the FUV-based SFR (e.g. Sullivan et al. 2004). However, this is not necessarily the case in low-luminosity 
dIrr galaxies that may have more fiuctuating SFHs. 

Figure [14] shows how well the GALEX FUV passband can recover the SFR averaged over the past 0.1 Gyr from 
our mock SFHs. The formula for transforming FUV to SFRpuv used here was from Hunter et al. (2010), where the 
Kennicutt (1998) formula was modified to account for the low metallicity of dwarf galaxies. Here we further divided 
their proportionality constant by 1.59 to account for the difference between the initial mass function (IMF) assumed 
here (Chabrier 2003) and the Salpeter (1955) IMF that they used, figure [15] justifies the practice of using FUV as an 
unbiased estimator of the recent SFR, with an uncertainty less than a factor of 2 incurred by different SFHs. One 
should keep in mind that here we assumed low dust extinction, which is reasonable for dwarf galaxies in general. 

Bias introduced by an incomplete SFH library 

The incompleteness of the SFH library used in the SED modeling could significantly affect or bias the estimation of 
physical parameters (e.g. Kauffmann et al. 2003). To demonstrate this, we create another SFH library by assuming the 
real SFH can be approximated as an underlying component that varies exponentially with time plus a single random 
burst of finite length. Here we denote this SFH library as 'E-l-B'. E-l-B is now the most commonly used SFH library 
in the literature. The exponential component is described with two parameters: the SF timescale r and the age. The 
burst component is described with three parameters: the age, strength, and length of the burst. The Chariot & Fall 
(2000) extinction recipe was adopted here. All the parameters, including the extinction, are allowed to vary uniformly 
within physically reasonable ranges when creating the library. We model the same SEDs as related to the mock SFHs 
generated above. We obtained the most probable physical parameters, including M^, SFRi, SFRq.i and SFR0.03, with 
the same method as described in § 14.21 

Figures [15] and [16] show the modeling results of using the E-f B library. Compared to the modeling results from 
our standard 6-C superposition method, the estimate of the relevant parameters exhibits a larger uncertainty and a 
considerable bias. For instance, the stellar mass M* has a tendency to be slightly underestimated, and the recent SFR, 
i.e. SFRq.i and SFRi, tend to be overestimated. The problem of using the E-l-B library in modeling can be more clearly 
seen in Figure [TBI The E-l-B library tends to overestimate the SFR averaged over recent times, e.g. SFRq.i, SFRi. As 
is expected for the more or less continuous nature of the SFHs in the E-l-B library, there is a very strong bias towards 
almost constant recent SF. The above bias could be ascribed to the fact that, while the SF averaged over a timescale 
of ~ Gyr may be relevant to the overall evolution histories of the galaxies, the SF over shorter timescales (i.e. 0.1 Gyr) 
is not necessarily as important. 



Outside-in Shrinking of the Star-forming Disk of Dwarf Irregular Galaxies 



15 



TABLE 1 
Galaxy Sample 
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Note. — (1) Galaxy names adopted in this work. (2) The other commonly used names in the literature. (3) Distance from Hunter ct al. , in prep 
and references therein. (4) Minor-to-major axis ratio measured on the V-band images. (5) The foreground reddening from Schlcgel ct al. (1998). 
(6) B-band absolute magnitude. (7) Integrated star formation rate (SFR,) normalized to the area of the galaxy within one V-band disk scale length. 
The SFR is derived from Ha. (8) The baryonic (stellar plus atomic gas) mass. The stellar mass is derived from our SED modeling in this work. 
The atomic gas (1.34xMffj) mass is collected from single dish observations in the literature (see Hunter & Elmegrcen (2004) for the references). 
(9) The atomic gas-to-stellar mass ratio. (10) Tidal index from Karachentsev et al. (2004). quantifies the collective gravitational disturbance 
from neighboring galaxies. The larger the value of is, the stronger gravitational disturbance exerted by the neighboring galax:ies. Galax:ies with 
zero or negative values of could be considered as isolated objects. 

^ The average of the E(B-V) derived by Burstein & Heiles (1982) (0.508) and by Schlegel et al. (1998) (0.700). See Johnson et al. in prep. 
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TABLE 2 
Multi-band Disk Scale Length 
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.31±0. 


,00 


0, 


.38±0, 


.00 


0. 


,46±0. 


,00 


0, 


.63±( 


NGC 6822 


1.29±0.00 


0.74±0.00 


0.63±0.00 


0.59±0.00 


0.65±0.00 


0.65±0.00 
























DDO 53 . . . 


0.20±0.00 


0.24±0.00 


0.31±0.00 


0.43±0.00 


0.54±0.00 


0.35±0.00 
























DDO 154 . 


0.60±0.00 


0.56±0.00 


0.52±0.00 


0.51±0.00 


0.52±0.00 


0.53±0.00 
























DDO 87 . . . 


0.56±0.01 


0.71±0.01 


1.94±0.01 


1.84±0.01 


1.73±0.00 


1.75±0.01 


0. 


.87±0. 


,03 


1, 


.40±0, 


.02 


0, 


,81±0, 


,01 


1, 


.03±l 


DDO 52 . . . 


4.30±0.03 


3.07±0.03 


1.39±0.00 


1.45±0.00 


1.33±0.00 


2.09±0.00 


0, 


,85±0, 


,03 


0, 


,89±0, 


,01 


0, 


,72±0, 


,07 


0, 


,93±l 


DDO 168 . 


-10.15±0.36 


-15.31±0.17 


174.25±19.25 


-58.18±0.89 


82.50±0.98 


-17.72±0.34 


0, 


,60±0, 


,00 


0, 


,60±0, 


,00 


0, 


,66±0, 


,00 


0, 


,69±( 


NGC 1569 


0.22±0.00 


0.23±0.00 


0.22±0.00 


0.22±0.00 


0.24±0.00 


0.31±0.00 


0, 


,05±0, 


,17 


0, 


,16±0, 


,58 


0, 


,60±0, 


,12 


1, 


,99±i 


NGC 3738 


0.21±0.01 


0.24±0.00 


0.36±0.00 


0.42±0.00 


0.45±0.00 


0.45±0.00 
























DDO 50 . . . 


1.78±0.01 


1.67±0.00 


1.37±0.00 


1.40±0.00 


1.27±0.00 


0.99±0.00 


0, 


,80±0, 


,01 


0, 


,75±0, 


,01 


0, 


,76±0, 


,01 


0, 


,82±i 


NGC 2366 


1.49±0.01 


1.39±0.00 


1.50±0.00 


1.58±0.00 


1.58±0.00 


1.83±0.00 


1, 


,07±0, 


,03 


1, 


,06±0, 


,01 


1, 


,08±0, 


,00 


1, 


,ll±i 


NGC 4214 


0.49±0.00 


0.46±0.00 


0.50±0.00 


0.55±0.00 


0.57±0.00 


0.54±0.00 


0, 


,85±0, 


,00 


0, 


.79±0, 


.00 


0, 


,71±0, 


,00 


0, 


,69±i 


NGC 1156 


0.68±0.00 


0.78±0.00 


0.88±0.00 


0.94±0.00 


0.96±0.00 


1.04±0.00 


0. 


.42±0. 


,02 


0, 


.52±0, 


.02 


0. 


,62±0. 


,00 


0, 


.72±( 



Note. — (1) Galaxy names. (2-7) Inner disk scale lengths measured on the surfaee brightness profiles of FUV, NUV, U, B, V., and 3.6 /J.m, 
respectively. The typical uncertainties of the scale length are smaller than 0.01 kpc. (8-13) Disk scale lengths measured in the outer part of the 
disks, if an obvious break is present in the surface brightness profiles. 
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TABLE 3 

SED-FITTING RESULTS 



Galeixy 


M, 


(SFR)l3.7Gyr 


R 


*^ 






log(E*, center) 


C31,* 


fco.lGyr 


blGyr 








Inner 


Outer 
















[10'^ Mol 


[10^^ Moyr-i] 

L ^ J 


(kpc) 


(kpc) 


(kpc) 


(log(Mopc-2)) 


(log(M.:.pc-2)) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


LGS 3 . . . . 


0.027 


0.038 


0.49 


0.09 


0.26 


-0.32 


-0.20 


2.04 


0.10±0.04 


0.17±0.06 


DDO 210 . 


0.068 


0.094 


0.44 


0.19 


0.29 


-0.21 


-0.05 


2.63 


1.72±0.64 


0.43±0.14 


DDO 69 . . . 


0.102 


0.137 


1.18 


0.23 


0.25 


0.03 


0.04 


2.36 


3.62±0.56 


1.55±0.30 


DDO 155 . 


0.297 


0.404 


0.22 








0.87 


3.06 


3.87±0.98 


0.84±0.21 


DDO 216 . 


1.521 


2.140 


0.89 


0.31 


1.12 


0.06 


0.44 


2.25 


0.06±0.02 


0.09±0.02 


MSldwA .. 


0.138 


0.182 


-14.28 


0.30 


0.41 


-0.12 


-0.16 


2.02 


7.81±1.12 


1.22±0.18 


Mrk 178 . . 


1.110 


1.531 


0.27 


0.87 


1.58 


-1.14 


0.69 


2.78 


2.04±0.52 


0.52±0.12 


DDO 187 . 


0.345 


0.451 


0.53 


0.23 


0.30 


0.48 


0.71 


2.51 


0.55±0.20 


2.88±0.62 


UGC 8508 


0.764 


1.061 


0.57 


0.27 


0.46 


0.50 


0.77 


2.49 


1.05±0.29 


0.56±0.16 


NGC 4163 


2.708 


3.807 


0.43 


0.29 


0.61 


0.83 


1.40 


2.62 


0.51±0.20 


0.32±0.10 


CVnIdwA . 


0.410 


0.532 


0.68 








0.22 


2.53 


1.57±0.41 


1.43±0.45 


DDO 70 . . . 


1.960 


2.668 


0.50 


1.59 


0.77 


0.14 


0.76 


2.77 


1.03±0.22 


1.44±0.25 


IC 1613 ... 


2.920 


4.030 


0.74 


0.39 


1.98 


-0.20 


0.96 


2.64 


0.81±0.25 


0.83±0.21 


DDO 101 . 


6.543 


9.325 


0.68 


0.41 


1.54 


0.39 


1.37 


2.52 


0.08±0.02 


0.12±0.05 


WLM 


1.629 


2.268 


1.24 


0.57 


1.34 


-0.08 


0.15 


2.31 


0.43±0.11 


0.27±0.07 


DD075 ... 


0.784 


1.023 


0.64 


0.27 


0.73 


0.31 


0.71 


2.03 


9.17±1.49 


2.21±0.31 


VIIZW 403 


1.777 


2.449 


0.96 


0.55 


1.14 


0.17 


0.67 


2.45 


3.51±0.89 


1.10±0.26 


Haro 29 . . . 


1.443 


1.921 


0.20 


1.98 


1.12 


-0.81 


1.59 


5.29 


6.88±2.02 


1.70±0.39 


DDO 133 . 


3.042 


4.214 


0.91 


0.57 


2.26 


-0.28 


0.77 


2.54 


0.72±0.14 


0.74±0.14 


DDO 126 . 


1.605 


2.132 


0.82 








0.37 


2.58 


1.42±0.33 


1.71±0.36 


DDO 165 . 


3.453 


4.631 


1.82 


0.92 


1.34 


0.22 


0.45 


2.30 


3.84±0.83 


1.34±0.25 


DDO 63 . . . 


3.337 


4.421 


1.61 


1.00 


1.57 


0.14 


0.56 


2.29 


2.25±0.52 


2.27±0.42 


NGC 6822 


7.630 


10.451 


1.47 


0.36 


0.38 


1.59 


1.66 


2.35 


0.51±0.14 


0.63±0.14 


DDO 53 . . . 


0.970 


1.337 


20.87 


0.79 


0.96 


-0.06 


-0.15 


2.10 


1.13±0.18 


0.69±0.14 


DDO 154 . 


0.835 


1.131 


1.00 


0.49 


0.68 


0.10 


0.28 


2.47 


3.82±0.83 


0.91±0.17 


DDO 87 . . . 


3.306 


4.589 


2.27 


1.06 


1.21 


0.15 


0.32 


2.69 


0.64±0.15 


0.53±0.13 


DDO 52 . . . 


5.273 


7.296 


1.05 


0.52 


3.15 


-0.41 


0.85 


2.68 


0.46±0.17 


0.86±0.25 


DDO 168 . 


5.865 


7.953 


32.04 


0.99 


0.73 


0.65 


0.63 


2.64 


1.43±0.29 


0.79±0.17 


NGC 1569 


36.016 


49.259 


0.43 


0.87 


2.25 


0.13 


2.35 


3.13 


6.19±1.36 


1.02±0.18 


NGC 3738 


46.571 


65.325 


0.68 


0.61 


1.59 


1.16 


2.20 


2.95 


0.62±0.22 


0.62±0.17 


DDO 50 . . . 


10.725 


14.416 


1.08 


0.71 


2.00 


0.42 


1.02 


2.45 


3.88±0.51 


1.47±0.21 


NGC 2366 


6.954 


9.186 


3.69 


1.11 


2.66 


-0.14 


-0.18 


2.70 


2.64±0.44 


2.00±0.30 


NGC 4214 


47.999 


65.390 


0.82 


0.60 


2.73 


0.56 


2.01 


3.09 


0.92±0.15 


1.26±0.22 


NGC 1156 


121.759 


164.914 


0.97 


0.64 


3.92 


0.55 


2.30 


2.64 


1.73±0.40 


0.85±0.22 



Note. — (1) Galaxy names. (2) Total stellar mass derived by applying the curve-ol'-growth method to the stellar mass profiles from our SED 
modeling. (3) SFR averaged over the Hubble time (13.7 Gyr). (4-5) Disk scale length measured from the surface stellar mass density profiles. The 
scale lengths of the inner and the outer parts are listed, respectively. If there is no obvious break in the profile, the related scale length is listed 
in column 4. (6) Radius where the break of the surface stellar mass density profile happens. (7) Surface stellar mass density at the broken radius. 
^*,(>reafe ^as been inclination-corrected under the assumption that the intrinsic (6/a)o = 0.3. (8) Inclination-corrected central surface stellar mass 
density, extrapolated from the inner regions around the center. (9) Concentration index (see the text for details). (10) Ratio of the SFR averaged 
over the past 0.1 Gyr to the SFR averaged over the Hubble time. (11) Ratio of the SFR averaged over the past 1 Gyr to the SFR averaged over the 
Hubble time. 



18 



Zhang et al. 



TABLE 4 

Slopes op the SFH Variations with Radius 



^ , A(log(SsFa„ A(log(SsFRo i/SQ) A(log(SsFRi /S,)) A(log(SsFRi 

^^^^y A(Rlkpcl) A(R/Rg) A(R|kpcl) A(R/R^) 

Inner Uuter Inner Uuter inner Uuter inner Outer 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



Galaxies with Mj,arj/ < 10® M© 



LGS 3 . . . . 


-0.10±0.16 








-0.02±0.04 








-1.25±0.19 




-0.29±0.04 








DDO 210 . 


-2.59±0.59 








-0.43±0.10 








-0.85±0.22 




-0.14±0.04 








DDO 69 . . . 


-0.50±0.06 


-2 


29±0 


28 


-0.09±0.01 


-0 


43±0 


05 


1.01±0.15 


-0.95±0.05 


0.19±0.03 


-0 


18±0 


01 


DDO 155 . 


-1.41±0.06 


-2 


53±0 


24 


-0.21±0.01 


-0 


37±0 


03 


-0.65±0.15 


-1.19±0.26 


-0.10±0.02 


-0 


18±0 


04 


DDO 216 . 


-0.44±0.04 








-0.24±0.02 








-0.52±0.04 




-0.28±0.02 








MSldwA . . 


0.58±0.03 


-0 


75±0 


02 


0.15±0.01 


-0 


20±0 


01 


0.56±0.03 


-0.65±0.01 


0.14±0.01 


-0 


17±0 


01 


Mrk 178 . . 


-2.53±0.12 


-1 


86±0 


02 


-0.66±0.03 


-0 


48±0 


01 


-1.28±0.09 


-0.35±0.02 


-0.33±0.02 


-0 


09±0 


01 


DDO 187 . 


-1.01±0.08 


-2 


09±0 


05 


-0.18±0.01 


-0 


37±0 


01 


-0.64±0.01 


-1.42±0.06 


-O.lliO.Ol 


-0 


26±0 


01 


UGC 8508 


-0.58±0.02 


-2 


56±0 


07 


-0.15±0.01 


-0 


68±0 


02 


0.07±0.02 


-0.72±0.02 


0.02±0.01 


-0 


19±0 


01 


NGC 4163 


-0.74±0.02 


-0 


06±0 


02 


-0.74±0.02 


-0 


06±0 


02 


-0.55±0.01 


0.00±0.04 


-0.55±0.01 





00±0 


04 


CVnIdwA . 


0.24±0.07 


-2 


50±0 


15 


0.14±0.04 


-1 


42±0 


09 


0.07±0.03 


-1.40±0.18 


0.04±0.02 


-0 


79±0 


10 


DDO 70 . . . 


-1.58±0.03 


-0 


82±0 


11 


-0.75±0.01 


-0 


39±0 


05 


0.00±0.05 


-1.42±0.11 


0.00±0.02 


-0 


68±0 


05 


IC 1613 . . . 


0.11±0.05 


-1 


22±0 


03 


0.06±0.03 


-0 


72±0 


02 


0.04±0.01 


-0.67±0.02 


0.02±0.01 


-0 


39±0 


01 


DDO 101 . 


0.52±0.01 


-0 


24±0 


07 


0.48±0.01 


-0 


23±0 


06 


0.68±0.03 


-0.18±0.08 


0.63±0.03 


-0 


16±0 


07 


WLM 


-0.46±0.02 


-0 


90±0 


04 


-0.26±0.01 


-0 


52±0 


02 


-0.26±0.01 


-0.39±0.06 


-0.15±0.01 


-0 


23±0 


03 


Averages . . 


0.01±0.72 


-1 


26±0 


68 


-0.10±0.31 


-0 


37±0 


17 


-0.29±0.65 


-0.71±0.35 


-0.14±0.16 


-0 


19±0 


07 



Galaxies with Mj,„^j, > 10® Mq 



DDO 75 . . . 


0.76±0.03 


-2.01±0.06 


0.17±0.01 


-0.45±0.01 


0.82±0.03 


-0.55±0.05 


0.18±0.01 


-0 


12±0 


01 


VIIZW 403 


-1.42±0.06 


-0.95±0.03 


-0.75±0.03 


-0.50±0.02 


-0.52±0.05 


-0.71±0.06 


-0.27±0.02 


-0 


37±0 


03 


Haro 29 . . . 


-0.60±0.09 


-1.09±0.05 


-0.18±0.02 


-0.32±0.01 


0.63±0.06 


-0.13±0.03 


0.18±0.02 


-0 


04±0 


01 


DDO 133 . 


0.46±0.01 


-0.23±0.05 


0.52±0.01 


-0.26±0.06 


0.15±0.01 


0.02±0.02 


0.17±0.01 





02±0 


02 


DDO 126 . 


-0.23±0.02 


-0.02±0.02 


-0.20±0.01 


-0.02±0.02 


0.42±0.03 


-0.02±0.01 


0.36±0.03 


-0 


01±0 


01 


DDO 165 . 


0.05±0.01 


-1.00±0.02 


0.12±0.01 


-2.27±0.04 


0.12±0.01 


-0.16±0.01 


0.26±0.02 


-0 


35±0 


02 


DDO 63 . . . 


-0.07±0.02 


-1.02±0.03 


-0.03±0.01 


-0.51±0.02 


0.19±0.03 


-0.39±0.02 


0.09±0.02 


-0 


19±0 


01 


NGC 6822 


-0.07±0.02 


1.75±0.14 


-0.04±0.01 


1.00±0.08 


-0.67±0.02 


0.44±0.13 


-0.38±0.01 





25±0 


07 


DDO 53 . . . 


-0.48±0.02 


-1.56±0.05 


-0.35±0.01 


-1.13±0.04 


-0.54±0.02 


-1.32±0.03 


-0.39±0.02 


-0 


95±0 


02 


DDO 154 . 


0.03d=0.03 


-0.45±0.02 


0.02±0.02 


-0.27±0.01 


-0.02±0.01 


0.03d=0.02 


-O.OliO.Ol 





02±0 


01 


DDO 87 . . . 


0.20±0.03 


-0.04±0.01 


0.29±0.04 


-0.06±0.02 


-0.05±0.01 


0.12±0.01 


-0.07±0.01 





17±0 


01 


DDO 52 . . . 


0.19±0.02 


0.13±0.01 


0.25±0.02 


0.17±0.02 


0.31±0.01 


-O.lOiO.Ol 


0.40±0.02 


-0 


14±0 


01 


DDO 168 . 


-0.06±0.01 


-0.88±0.02 


-0.05±0.01 


-0.73±0.02 


0.01±0.02 


-0.40±0.02 


0.01±0.02 


-0 


33±0 


02 


NGG 1569 


-0.64±0.03 




-0.25±0.01 




-0.46±0.03 




-0.18±0.01 








NGC 3738 


-1.31±0.04 


-0.27±0.01 


-1.31±0.04 


-0.27±0.01 


-1.22±0.02 


-0.25±0.02 


-1.22±0.02 


-0 


25±0 


02 


DDO 50 . . . 


0.35±0.01 


-0.31±0.02 


0.38±0.01 


-0.34±0.03 


0.08±0.01 


O.OOitO.Ol 


0.09±0.01 





00±0 


01 


NGC 2366 


0.04±0.01 


-0.07±0.01 


0.06±0.01 


-0.10±0.02 


-0.14±0.01 


-0.09±0.01 


-0.19±0.02 


-0 


12±0 


02 


NGC 4214 


0.15±0.01 




O.llitO.Ol 




-O.OlitO.Ol 




O.OOitO.Ol 








NGC 1156 


-0.30±0.01 


-0.41±0.01 


-0.24±0.01 


-0.33±0.01 


-0.05±0.01 


-0.25±0.02 


-0.04±0.01 


-0 


20±0 


02 


Averages . . 


0.12±0.29 


-0.59±0.44 


0.08±0.25 


-0.44±0.65 


0.04±0.21 


-0.08±0.20 


0.02±0.20 


-0 


09±0 


16 



Note. — (1) Galaxy names. (2-5) Slopes of the radial variations of A(log(I]sFRo.i/S*))/A(R[kpc]) and A(log(i:sFRo.i/2*))/A(R/R^), for the 
inner disks and the outer disks, respectively. (6-9) Slopes of the radial variations of A(log(EsFRj^ /Xl*))/A(R[kpc]) and A(log(SgFRi/S^,))/A(R/R;p), 
for the inner disks and the outer disks, respectively. 
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Fig. 1. — First column: Azimuthally-averaged surface brightness profiles of tfie Ha, FUV, B and 3.6 fivn images. Tfie NUV or U-hand 
profiles are plotted instead of FUV if there are no FUV observations. Second column: Radial variations of the SFR averaged over different 
timescales obtained from our multi-band SED modeling. The averaging timescales are the most recent 0.03 Gyr, the most recent 0.1 Gyr, 
the past 1 Gyr, and the Hubble time. Third column: Inclination-corrected stellar mass surface density profiles from the SED modeling. 
The fitted exponential profiles are overplotted as solid lines. The thin red line marks the Holmberg radius, and the thick black vertical line 
marks the V-band scale length. [A color version of this figure is available in the electronic journal.] 
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Fig. 1- Continued 
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RB /t-) 3.6ij/m 
D /^D 

Fig. 2. — Ratios of disk scale lengths /R^jj^^" plotted against 7?^//?^^^"° for our sample of galaxies. R™ , R^, and R^j^^" denote 
the disk scale lengths of the FUV, B, and 3.6 fiia passbands, respectively. The disk scale length measured at the outer disk is plotted for 
galaxies with broken surface brightness profiles. The galaxies with baryonic mass larger and smaller than 10® Mq are denoted as black open 
circles and red open triangles, respectively. The (solid) line of equality is shown to guide the eye, and galaxies below the line have larger 
than i?^. [See the electronic journal for a color version of this figure.] 




Fig. 3. — The disk scale length ratios of the outer- to-inner disks in the FUV plotted against the disk scale length ratios of the outer-to-inner 
disks in the 3.6 /xm. The solid line is the relationship of equality. The dotted lines mark values for non-broken profiles. 
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Fig. 4. — Fractions of the total stellar mass formed during the past 1 Gyr fioyr derived from our SED-modeling are compared to those 
determined from stellar CMD-based analysis. The galaxies with CMD analysis from Weisz et al. (2011) and from Orban et al. (2008) are 
denoted as filled circles and filled triangles, respectively. The {solid) line of equality is plotted to guide the eye. The dotted lines mark the 
expected value for a constant SFH over the Hubble time. 
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Fig. 5. — Ratios of globally integrated SFR0.1/SFR13.7 vs. SFR1/SFR13.7 for the galaxy sample. The {solid) line of equality is plotted to 
guide the eye. The galaxies are divided into two groups according to total baryonic mass. Those with masses less than IO^Mq are denoted 
with red triangles, and those with masses greater than this are plotted as black circles. [See the electronic journal for a color version of this 
figure.] 
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Fig. 6. — log(SsFRo.i/S*) is plotted as a function of radius, which is normahzcd by the K-band disk scale length. The galaxies are plotted 
in order of total baryonic mass from the upper left to the lower right, and within each panel galaxy names are listed in order of increasing 
baryonic mass. The dashed line in each panel marks a constant SFH over a Hubble time. 
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Fig. 7. — log{SsFRi /S*) plotted as a function of radius, which is normalized by the F-band disk scale length. The galaxies are plotted in 
order of total baryonic mass from the upper left to the lower right, and within each panel galaxy names are listed in order of increasing 
baryonic mass. The dashed line in each panel marks a constant SFH over a Hubble time. 
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Fig. 8. — Fitted slope of the radial variations of log(SsFRo plotted against the slope of the radial variations of log(SsFRi /S*). The 

left panel shows the results for the inner disks, and the right panel shows the results for the outer disks. The galaxies with baryonic mass 
larger and smaller than 10* Mq are denoted as black open circles and red open triangles, respectively. The large diamonds represent the 
averages among the whole sample galaxies. The (solid) line of equality is plotted to guide the eye. The dashed lines mark a fiat gradient 
slope. Here the slope fitting was done on the physical (kpc) scales. [See the electronic journal for a color version of this figure.] 
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Fig. 9. — Fitted slope of the radial variations of log(SsFRo plotted against the slope of the radial variations of log(SsFRi /S*). The 

left panel shows the results for the inner disks, and the right panel shows the results for the outer disks. The galaxies with baryonic mass 
larger and smaller than 10* Mq are denoted as black open circles and red open triangles, respectively. The large diamonds represent the 
averages among the whole sample galaxies. The (solid) line of equality is plotted to guide the eye. The dashed lines mark a fiat gradient 
slope. Here the slope fitting was done on the radius normalized to the K-band disk scale length. [See the electronic journal for a color 
version of this figure.] 
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Fig. 10. — Fitted slope ol the radial variations of log(SsFRo.i /S*) and log(SsFRi /S*) in the outer disks plotted as a function of the tidal 
index 0. The black circles and the red triangles denote galaxies with baryonic mass larger and smaller than 10* Mq, respectively. The open 
and the filled symbols denote the slope of log(SsFRo.i /S*) and the slope of log(EsFRi /S*), respectively. The dashed lines mark the upper 
and the lower envelopes of the distribution of most data points in the plot. [See the electronic journal for a color version of this figure.] 



Outside-in Shrinking of the Star-forming Disk of Dwarf Irregular Galaxies 



35 




-2 2 4 

log(S. [Mo pc-*]) 



Fig. 11. — SsFRo 1 plotted against for the galaxy sample. The data points are from the azimuthal averages of the relevant surface 
density profiles, so each galaxy is represented by majiy points. The solid line in the plot indicates the expected relationship between SsfRo.i 
and S» for a constant SFR over a Hubble time. 
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Fig. 12. — Histograms of the ratios of the parameters (M*, SFRi, SFRo.i and SFR0.03) from the input mock SFHs to the best estimates 
derived from the 6-component (6-C) modeling. The abscissa in eaeh panel is the logarithm of the relevant ratio (indicated in each panel) of 
the input and of the 6-C modeling results. The ordinate is the per cent of the data points in each bin. The dotted line in each panel marks 
the relationship of equality. The red curve in each panel is the gaussian fit to the relevant histogram, cr and fi indicate the best fitting 
parameters defining the gaussian curve. 
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Input SFHs 



Fig. 13. — Comparison between the ratios of the SFR averaged over different timescales (indicated in each panel) from the input SFHs 
(the abscissa) and the relevant quantities from the 6-C modeling (the ordinate). The red solid line in each panel marks the relationship of 
equality. 
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Fig. 14. — Histogram of the ratio of SFRo.i from the input mock SFHs to SFRfuv estimated with the modified formula of Hunter et al. 
(2010). The ordinate is the per cent of the data points in each bin. The red curve is a single gaussian fitting to the histogram, cr and fi 
indicate the best fitting parameters defining the gaussian curve. The dotted line marks the relationship of equality. 
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Fig. 15. — Histograms of the ratios of the parameters (M*, SFRi, SFRo.i and SFR0.03) from the input mock SFHs to the best estimates 
from the exponential plus burst (E+B) modeling. The abscissa in each panel is the logarithm of the relevant ratio (indicated in each panel) 
of the input and of the E+B modeling results. The ordinate is the per cent of the data points in each bin. The red curve in each panel is 
the gaussian fitting to the relevant histogram, cr and fj. indicate the best fitting parameters defining the gaussian curve. The dotted line in 
each panel marks the relationship of equality. 
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Input SFHs 



Fig. 16. — Comparison between the ratios of the SFR averaged over different timescales (indicated in each panel) from the input SFHs 
(the abscissa) and the relevant quantities from the E+B modeling (the ordinate). The red solid line in each panel marks the relationship 
of equality. 



